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We discuss the thermodynamics of the O(N) model across the corresponding phase transition 
using the two-loop ^-derivable approximation of the effective potential and compare our results to 
those obtained in the literature within the Hartree-Fock approximation, fn particular, we find that 
in the chiral limit the transition is of the second order, whereas it was found to be of the first order in 
the Hartree-Fock case. These features are manifest at the level of the thermodynamical observables. 
We also compute the thermal sigma and pion masses from the curvature of the effective potential, 
fn the chiral limit, this guarantees that Goldstone's theorem is obeyed in the broken phase. A 
realistic parametrization of the model in the N = 4 case, based on the vacuum values of the 
curvature masses, shows that a sigma mass of around 450 MeV can be obtained. The equations 
are renormalized after extending our previous results for the N = 1 case by means of the general 
procedure described in Ref. [ij. When restricted to the Hartree-Fock approximation, our approach 
reveals that certain problems raised in the literature concerning the renormalization are completely 
lifted. Finally, we introduce a new type of <f>-derivable approximation in which the gap equation is 
not solved at the same level of accuracy as the accuracy at which the potential is computed. We 
discuss the consistency and applicability of these types of "hybrid" approximations and illustrate 
them in the two-loop case by showing that the corresponding effective potential is renormalizable 
and that the transition remains of the second order. 
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I. INTRODUCTION 

It is a well-known fact that the ^-derivable approxi- 
mation scheme, also called in the literature two-particle 
irreducible (2PI) or Cornwall- Jackiw-Tomboulis (CJT) 
formalism, gives a first order phase transition when ap- 
plied to the O(N) model in its lowest, Hartree-Fock ap- 
proximation level. Other resummation methods, such 
as the 1/N expansion give a second order phase transi- 
tion already at leading order @] , in accordance with gen- 
eral expectations and with the result of the functional 
renormalization group approach 0, Qj. It was argued 
0, ||| that close to the transition temperature the contri- 
bution of higher loops may become important and that 
already the inclusion of the setting-sun diagram in the 
^-derivable functional will render the phase transition of 
the second order type. 

As a continuation of our previous investigation done for 
the one-component scalar model, where we found that 
the change of order indeed happens within a full two- 
loop treatment of the effective action, we turn now to 
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the physically more interesting 0(N) model. For N — 
4 this model can be regarded as a low energy effective 
model of two flavor QCD because the global 5C/(2)l x 
SU (2)r symmetry of the latter is isomorphic with 0(4). 
Since the 0(4) model contains both the longitudinal and 
transverse excitations of the chiral order parameter, it is 
widely used in the phenomenology of low energy mesons 
for the qualitative description of medium induced effects, 
especially around the phase transition. 

We would like to understand where exactly the contri- 
bution coming from the setting-sun diagram is essential 
to obtain the right order of the phase transition. There- 
fore, in addition to the two-loop approximation, we con- 
sider an approximation where the effective action is com- 
puted at two-loop order, but is evaluated for propaga- 
tors computed from the Hartree-Fock approximation. Al- 
though this hybrid type of approximation might present 
certain inconsistencies, as it does not obey the conditions 
identified by Baym in Ref. Q , it is convenient in practice 
because its numerical treatment is much easier. More- 
over, we will see that, for not too low temperatures where 
the hybrid approximation does not seem to show incon- 
sistencies, its results are pretty close to those obtained 
from the two-loop approximation which is numerically 
more time and memory demanding. In particular, both 
in the two-loop and in the hybrid approximations, the 



2 



transition is found to be of the second order. 

As far as meson phenomenology is concerned, we will 
be particularly interested in the value of the sigma mass 
which can be obtained within a realistic paramctrization 
of the model. Indeed, one of the difficulties when apply- 
ing the 0(4) model to meson phenomenology is to obtain 
high enough values of the sigma mass, while maintaining 
the interpretation as an effective model where the cutoff 
A, a mere separation scale between the modes of interest 
(p <C A) and those which arc integrated out (p > A), does 
not play the role of a parameter. This is usually rendered 
difficult by the fact that the model possesses a Landau 
pole in the ultraviolet and, if the latter is too close to 
the physical scales, the rcnormalization procedure is not 
enough to ensure the insensitivity of the results to cut- 
off values below the Landau pole. We will see that, in 
the two-loop and hybrid approximations, one can obtain 
reasonable values of the sigma mass with a Landau scale 
almost one order of magnitude higher. This, combined 
with the fact that the Landau pole does not show up in 
the renormalizcd quantities defined within the two-loop 
or hybrid approximations, allows us to meet the above 
mentioned requirements. 

In principle the insensitivity to the cutoff scale A is en- 
sured automatically by the renormalization group since, 
by following a line of constant physics, a change in A is 
carried over to the (bare) parameters of the Lagrangian, 
in such a way that the low energy physics is unaffected. 
Even if this picture persists order by order in pertur- 
bation theory, this is not necessarily so for approxima- 
tion schemes that go beyond it and certain amendments 



need to be made to the renormalization procedure, de- 
pending on the method used. Over the last few years 
a general method for renormalizing ^-derivable approxi- 
mations has been developed and we illustrate it here both 
in the two-loop and in the hybrid approximations to the 
2PI effective action. By revisiting the lower Hartree-Fock 
approximation, wc can also compare our renormalization 
procedure to other approaches followed in the literature. 
In particular, we show that certain inconsistencies dis- 
cussed in Ref. [|[ are completely lifted by our approach. 

In Sees. |TT] and IIII1 we define and renormalize the 
two approximations to be discussed in this work and 
compare our renormalization procedure to other ap- 
proaches. Section Hvl is devoted to some of the numerical 
tricks that we use to achieve high accuracy results in 
the two-loop approximation. Section [V] deals with the 
paramctrization of the model, with a special attention 
to the attainable values of sigma mass and gathers our 
results on the phase transition and thermodynamical 
observables. We also discuss there the dependence of 
the physical quantities on the renormalization scale and 
on the cutoff. We present our conclusions in Sec. IVII 

II. TWO-LOOP APPROXIMATION 

A. Relevant equations 

The 2PI effective action for the O(N) model is a func- 
tional of a one-point function (f> a (x) and a symmetric two- 
point function G a b(x,y) = Gba(y,x). In the imaginary- 
time formulation of field theory at a finite temperature 
T = 1/(3 and at two-loop order, it reads 
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(1) 



with f x = J Q dr J d 3 x, 4> 2 = a a , 4>G<j) = <p a G a b<pb, 
tr G = G aa and where a summation over repeated 
indices is implied. As discussed in Rcfs. [l|, Q and 
below, the need for two bare masses mo and m-2 and 
three families of bare couplings, labeled with the indices 
0, 2 and 4 respectively, reflects the fact that, given a 
truncation of the 2PI effective action, there are two 
possible ways to define the two-point function and three 
different ways to define the four-point function. As 
we discuss in Appendix \K\ the need for a doubling 
(represented by the superscripts A and B) of the bare 
couplings carrying an index or 2 has to do with the fact 
that two of these four-point functions do not obey the 



crossing symmetryQ Finally, following Ref. we have 
replaced the bare couplings in the highest loop diagrams 
of Eq. ([T]) by a coupling A* which will be identified later 
to the renormalizcd coupling at some renormalization 
scale This is because no renormalization comes from 
these vertices at this level of truncation. 



1 Equivalently, the corresponding terms in the 2PI effective action 
JTJ are independently invariant under O(N) transformations, see 
Ref. Q. 
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In what follows, we study the phase transition of the 
model by computing the effective potential 7(0). The 
latter is obtained after evaluating the functional (fTJ) at 
the stationary value of G which we denote G^0 with 4> 
a homogeneous field configuration: 



7(0) =-pyT[<t>,G 4 



(2) 



Some more explicit expressions of the effective potential 
will be given later. In the presence of a homogeneous 
field, the propagator G(x, y) depends on the difference 
x — y or, in Fourier space, on Q = (iw n ,q) where lu„ = 
27rn//3 is a bosonic Matsubara frequency. From parity 
and time-reversal symmetry, we have G ab (Q) = G ab (—Q) 
and thus G ab (Q) = G ba {Q) since G ab (Q) = G ba (-Q). 
Moreover, the O(N) invariance of T[<fi, G] upon simulta- 
neous rotation of <f> and G implies that G^ is covariant: 

G R J = R ac RbdGi d , VReO(N). (3) 

We recall in Appendix [C] that, together with the prop- 
erty G ab (Q) = G ba (Q), this implies the following spectral 
decomposition 



G ab ~ GhPah + G T Pj b , 



with 



and Pju 



(4) 



(5) 



the longitudinal and transverse projectors with respect 
to <j) an d where the functions Gl and Gt depend on cj> 
only through </> 2 . 

It is convenient to introduce momentum dependent 
longitudinal and transverse masses defined through the 

relation G l ,t(Q) = l/(Q 2 + M'l . T (<2))> such that thc Y 
include the corresponding self-energy and the tree-level 
mass. After some straightforward calculation starting 
from the stationarity condition = ST/SG\q, one shows 
using the two projectors in Eq. ([5]) that they obey the 
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following coupled gap equations! 



Ml(K)=m 2 + 
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67V 



-T[G L 



A 



((JV-l)A) 




67V 



-T[G4 



^'■[9B[G L }(K) + (N-1)B[G T }(K)} 



A 



(A+2B) 



67V 



(0) 



and 



a(^) ,((N-1)A+2B) 

M$(K) =m 2 + — — T[Gl] + ^ 

TV ' 67V 1 1 67V 



T[G T ] 



2 In order to alleviate the notations, the dependence of G and alike 
on various quantities such as <j>, T, . . . will be written explicitly 
only when needed. 

3 By using the substitutions Aq^~*~^ S / (a + f3) — > 12/e and A* — > 
12/e, we recover the equations derived in Ref. lid . There how- 
ever, the equations were neither renormalized nor solved. 
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\i A) -^B[G L ;G T }(K) 



(7) 



In order to save space, we find it appropriate to write 



A 0,2 = aA Q,2 + P A 0,2 i 



and denoting the sum integral by 

<- T ' d*q 

;wj = j 1 ■ 

IQ 

where q = \q\, we use the short-hand notations 

i-T 



r T 00 r d 3 n 

i /(Q)3r l/R ,M 



(8) 



(9) 



T[G] ee f G(Q) , (10) 

B[G 1 :G 2 ](K)= f G 1 (Q)G 2 (Q + K) , (II) 
Jo 



S[G 1 ;G 2 ;G 3 }= [ f G 1 (K)G 2 (Q)G 3 (Q + K) . 

JK JQ 

(12) 

For the last two of them, when all the arguments are 
equal to a given propagator G, we write more simply 
B[G]{K) and S[G\. Taking the difference of Eqs. © and 
(O, it is straightforward to check that, when = 0, the 
system of equations is compatible with a solution such 
that Ml = M| ee AfJ =0 with 



(NA+2B) 



67V 



-T[G, 



(13) 



and G 0=o (O) = 1/(Q 2 + M2 =o ). 

The nature of the transition will be discussed by moni- 
toring the nontrivial cxtrema 4> of the effective potential. 
They obey the equation 



= m + 
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67V' 
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67V 



-T\G^ 



A 



^(35[G l ] + (A^-1)5[G l ;Gt;Gt]), (14) 



187V 

which, due to the stationarity condition = ST/SG\q, 
originates only from the explicit field dependence of the 
functional ([T]). We note that the case TV = 1 is obtained 
after disregarding Eq. ([7]) and making the replacements 

\ { N 2 A+2B) = \ { A 2 +2B} -> 3A ,2 in Eqs. © and We 
shall use this recipe later in order to cross-check the 
expressions obtained for the bare parameters. 

We shall also need the curvature of the potential, which 
at d> — is found to be 



ANA+2B) N „ 

^ 1 N + 2 xls[G, =0 ] 



M 4>=o = 



67V 



-T[G0 = o] 



187V 2 



(15) 



More generally, one can define the curvature tensor at an 
arbitrary value of the field. Using as in Ref. [ll[ that the 
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effective potential depends on the field </> only through the 
0(iV)-invariant 2 , one writes 7(</>) = U((f> 2 ) and obtains 



5 2 7 (0) 



AU 



[2U'{4> 2 ) + 4^V(«/> 2 )]P L & + 2U'(cf> 2 )Pl 



(16) 

In this paper we shall call curvature masses the two eigen- 
modes appearing in the above equation, evaluated at the 
solution <j) of the field equation: 

Ml = 2U'(4> 2 ) +4<j) 2 U"(4> 2 ) and M| = 2[/'(0 2 ).(17) 

The field equation reads 



hit) 



2C/'(0 2 )0a =0, 



(18) 



from which it follows first, that M 2 = M 2 , in the symmet- 
ric phase (since <f> = 0) and second, that = in the 
broken phase (since 4> 7^ and thus U'(4> 2 ) = 0) in agree- 
ment with Goldstone's theorem. In contrast, there is no 
reason for the gap mass M| = M^,(K = 0) to vanish 
in the broken phase and we shall investigate quantita- 
tively how much the Goldstone's theorem is violated in 
this case. 

In the case of explicitly broken symmetry, when a 
term — h(j> = -~h a <j) a is added to the effective potential, 
what changes is the field equation^ which becomes 
S r y(4>)/S(j) a \^ = 2U'((f> 2 )4>a — h a = 0, so that we have 
M| = ||ft||/||<£|| and Ml still given by Eq. ([I7J). Without 
loss of generality, we can choose h = (\\h\\,0, . . . ,0) 
along the first coordinate axis. Note also that if 
we view 2U'((f> 2 ) = /(||0||) as a function of ||</>||, 
we can compute the longitudinal curvature mass as 
Ml = + ||0||/'(||</>||) horn a numerical derivative 

of the function /(||0||) which appears on the right-hand 
side of Eq. (HI) . 



The gap and field equations ©, (JT)) and (|14[) will be 
solved using the techniques developed in Ref. @ that 
we quickly summarize in Sec. IIVI Before we proceed to 
the numerical resolution of the model, we must however 
determine the values of the bare parameters in such a 
way that the sensitivity to the ultraviolet regulator is re- 
moved, or at least considerably reduced. The results that 
we shall present are valid for any regularization that can 
be defined nonperturbativcly. For definiteness however 
and in line with the numerical method that we use to 
solve the two-loop approximation, in the next section, 
we assume that 3D momenta of modulus larger than a 
given cutoff A are dropped. More details concerning the 
regularization procedure can be found in Ref. [9| . 



4 The stationarity condition that defines is not changed. It 
follows that Eq. J3]l and in turn Eq. {2) still hold. 



B. Renormalization 

As explained in Ref. [l[ and illustrated in Ref. the 
fact that the gap masses at zero momentum are different 
from the curvature masses requires the presence of two 
distinct bare masses m and m 2 . Those are fixed by 
means of the usual renormalization condition 



M 2 =0 ^(K = 0) 



(19) 



at some renormalization scale, here a temperature T*, 
supplemented by a consistency condition 



M l=0,T t 



M 2 (K = 0). 



(20) 



which restores the equality of the two masses at the renor 
malization point. Applying these conditions to Eqs. (|13| 
and (fl"5"1). we obtain 



A, 



(NA+2B) 




6N 



(21) 



and 



A 



(NA+2B) 



6N 



-TJGJ 



N ■ 



18N'< 



A 2 S*[Gy, (22) 



with G*(Q*) = 1/(Q 2 + m 2 ). The 7k- on any quantity 
means that it is computed at the temperature T*. For 
instance, Q+ means that the corresponding Matsubara 
frequencies involve the temperature T*. Similar consid- 
erations apply to the four-point function which admits 
three distinct definitions V ab , c d, V a b,cd(K), V abcd ; see Ap- 
pendix [XJ The first two do not obey the crossing symme- 
try and thus involve two independent components at <j> = 

0: Kt.7c° = Vj^lSab&cd + vj=l{3ac$bd + S ad Sbc) and sim- 
ilarly for V^ d {K). In contrast = %=o{S a b8 cd + 
5 ac Sbd + SadSbc) is crossing symmetric. The renormaliza- 
tion condition 



V 



(A) 
=0,T* 



3N 



(23) 



and the consistency conditions 

t> _ fM) _ fAB) 

V<j>=0,T* ~ V <t>=0,T is - %>=0,T„ 

vi% t {K = Q) = V^{K = Q) (24) 



fix all the bare couplings that we have introduced and 
restore the equality of the various four-point functions 
at the renormalization point, in particular, the crossing 
symmetry becomes manifest. We obtain the following 
expressions for the bare parameters: 



3N 3N 



A, 



(B) 



A 



B*[G*](0) 



(25) 



and 



3N 



3N 



X (NA+2B) (JV + 2)A* 2 



--B*[G*](0), (26) 



5 



for those coupling parameters labeled with 0, A 
A^f+SA^ with 



(B) 



(27) 



and 



A 2 / G5(Q*)AB*[G*](Q*), (28) 



187V 2 



as well as Af A+2B) = A<f A+2B) + 5X£ A+2B) with 



c\(NA+2B) _ A/ + 2 

and 



A 2 ^[G*](0) (29) 



A 



l 21 

(NA+2B) 




1 



A 2 / T ' Gl(Q+)AB4G*](Q+) , 



6N 2 



Q* 



(30) 



for those coupling parameters labeled with 2 and finally 



A 4 



-2A, 



1 (A 



(NA+2B)} 
21 J 



A 



(NA+2B) 



N 

{N + 2f (N + 6) 2 



2 1- 



(A 



(B, 
21 



67V 4 



54iV 3 



1 



G 2 (Q*) [AB*[G*](Q*)] : 



In the above expressions, AS*[G*](Q*) stands for the dif- 
ference of bubble sum integrals B ir [G ir ](Q ir ) — B*[G*](0). 
The reason for the splitting of the bare parameters X^' 3 ^ 

into "local" and "nonlocal" parts, xif' B ^ and SX^^ 
respectively, is explained in Refs. {12H14I ]: see also Ap- 
pendix [51 Applying the replacement rule discussed right 
after Eq. lfT4"|) , one recovers the N — 1 bare parameters 
of Ref. [0] . It is also simple to obtain the expressions for 
the bare parameters in the Hartree-Fock approximation. 

One has simply to set 5X^1^ — and to remove all those 

terms that involve A£>* [G*] . 

unchanged, while to 2 = too 
and 



Then Xq A ^ and Xq B ^ remain 



A! 



U) 



= A 



(A) A (B) 



= A 



(B) 



A 4 H = -2A* + X { a A+2B) 



(32) 



which gives A^ = — 2A+ + 3Ao when N = 1, in agreement 
with the result of Ref. [l5| . We have introduced a super- 
script 'H' for the value taken by A4 in the Hartree-Fock 
approximation for later convenience. 

Following similar steps as in Ref. @, it is possible to 
prove implicitly that the bare parameters given above 
renormalize the gap and field equations, as well as 
the effective potential (up to a temperature and field 
independent divergence for this latter quantity). By 
"implicit proof" , we mean that certain steps require 



some assumptions on the properties of a function, the 
spectral function, which is defined implicitly. Wc are not 
able to prove these properties but we can argue that they 
are plausible for they are true perturbatively and the 
resummation should only bring innocuous logarithmic 
corrections to them. We shall not reproduce this proof 
here and refer to Ref. [i| for further details. In the 
next section however, we illustrate some of the aspects 
of the proof which are specific to the O(N) model by 
using a simpler approximation where renormalization 
can be performed in an explicit way. This will be also 
the opportunity to revisit the renormalization of the 
Hartree-Fock approximation from our point of view 
and to compare to other results in the literature, in 
particular those of Ref. H . 



C. Landau pole 

Let us end this section by discussing the presence 
of a Landau pole in the O(N) model and how this 
affects the discussion of renormalization at the level of 
approximation considered in this work. 

First of all, at least one pole is present in the expres- 
sions for the bare parameters. Indeed, the equations 



and (|2"6")) determining the bare couplings Xi A ^ and Aq B ' ) 



(31) 

can be rewritten as 



1 



2A* 
~6N 



B*[G*](0) 



and 



A 



(A) 



A 



(AT + 2)A, 
6N 



B*[G*](0) 



(33) 



(34) 



where we have made the cutoff dependence of the bubble 
sum integral explicit. Since the latter grows logarithmi- 
cally with A, it follows that both Xq A ^ and Xq B ^ diverge 
before turning negative at some value of A, which signals 
an instability. The bare coupling Xq A ^ being the first to 
diverge since A^ > 0, it is natural to define the Landau 
scale A p from the equation: 



= 1- 



6N 



-AA Ap [GJ(0) 



(35) 



Above this scale, at least one of the bare couplings 
becomes negative and one might wonder whether the 
theory is stable. In contrast, below this scale, it is 
easily checked, using the fact that B*[G*](0) > 
and AB+{G±](Q+) < (this is proven for instance in 
Appendix B.3 of Ref. j9|), that all the bare couplings 
remain positive. To remain in the stability region, we 
shall thus consider values of A below A„. 
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In the case of the two-loop approximation considered 
here (and also in the hybrid approximation that we intro- 
duce in the next section or in the Hartrcc-Fock approxi- 
mation), the presence of a pole in the cutoff dependence 
of the bare couplings does not imply the appearance of a 
pole in the integrals that enter the physical observables. 
Choosing parameters such that the Landau scale is not 
too close to the physical scales!! the physical quantities 
are defined for any value of the cutoff A. This is because 
in the two-loop approximation (and also in the Hartree- 
Fock approximation or in the hybrid approximation con- 
sidered in the next section) the self-energy does not grow 
quadratically at large frequency /momentum and also be- 
cause these approximations do not involve vertex-type re- 
summations capable of generating a Landau pole in the 
physical quantities. It follows that one can discuss renor- 
malization as usual, in terms of divergent and convergent 
quantities as A — > oo and thus, even though we restrict to 
values of A below A p , the renormalization procedure en- 
sures that the results are already pretty insensitive to the 
cutoff in this range if the Landau scale is large enough. 
We have already studied these features in Refs. [1, 
and we shall also do it here briefly in Sec. [V] 

At higher orders of approximation, one expects a pole 
to appear in the physical observables too, at a finite value 
of the cutoff. This prevents discussing the renormaliza- 
tion in terms of divergent and convergent quantities as 
A — > oo. Still, if the Landau scale is large enough, these 
concepts survive in a somewhat generalized acceptation. 



In particular, quantities rcnormalized according to our 
scheme will still show a plateau behavior below the Lan- 
dau scale, from which one can extract results that are 
pretty insensitive to the cutoff. The discussion becomes 
more delicate as the Landau scale gets closer to the phys- 
ical scales. 



III. HYBRID APPROXIMATION 

We shall also consider another type of approximation 
where the gap equation is solved at a lower level 
of accuracy than that used to compute the effective 
potential. We name these approximations "hybrid" 
for they break to some extent the consistency of the 
(f>-dcrivablc formalism. In particular, because the po- 
tential is not evaluated at its stationary point, the field 
equation admits additional contributions of the form 
8T/5G\q 5G/Scf>. These types of approximations have 
been considered in earlier works as well; see [l6l - [l8j |. 
We note that these types of approximations do not 
obey Baym's conditions and might thus lead to certain 
inconsistencies in some region of the parameters. 



A. Definition and relevant equations 

To make things explicit, we consider the two-loop 2PI 
effective potential (below Cl = 1 and ct = N — 1): 



7 [^G L ,G T ]=iV7o(m.) + im^ 2 + -^^+ ]T | f [in Gr\Q) - In G^(Q) + (Q 2 + rog) G t {Q) - l] 

i=T,L Q 



x(A+2B) ,((N-1)A) ,((N-1) 2 A+2(N-1)B) 



T 2 [G L ] + ^— - T[G L }T[G T ] + -2 ^ r 2 [G! 



24 jV 12 N J 24iV 

AA+2B)^-\ n , . ((TV-l^Wr^ A 2 2 



N [X ( 2 A+AB) T[G L ] + Ar - 1)A> T[G T ]] - ^ [35[G L ] + (JV — 1)5[G L ; G T ; G T ]] , (36) 



but instead of evaluating it at its stationary point, de- 
fined by the solution of Eqs. (jSJ) and ([7]), we evaluate 
it at the stationary point of the Hartree-Fock effective 
potential. There are two main reasons to do this here. 
Since the Hartree-Fock gap equations are equations for a 
momentum independent self-energy, the possibility rises 
to draw some conclusions, including renormalization, an- 
alytically and also numerical calculations become faster, 



allowing for a thorough investigation of the model@ The 
momentum independence of the self-energy also allows us 
to conveniently work in dimensional rcgularization. We 
stress however that what follows can be redone equiva- 
lently using a cutoff rcgularization. We note finally that 
the need for 7 (m^) in the expression (|36[) stems from 
a proper regularization of the 2PI effective action, as 
discussed in Ref. §. In dimensional regularization, we 



If the Landau scale is too close to the other scales, we have seen 
in Ref. [3 that the gap equation might lose its solution if the 
cutoff is taken too large, implying that the physical observables 
are not denned for too large values of the cutoff. But this is not 
due to the appearance of a pole in the integrals contributing to 
these observables. 



6 We shall see that, once a physical parametrization of the model 
is performed, our results in the two-loop and hybrid approxima- 
tions will not differ much. 
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havofl 

7o (m*) 



1 



2 y (2?r) d - 1 
with d = 4 - 2e. 



+ 2Tln(l-e-^ /T ) , (37) 



Although it can be performed explicitly, see below, the 
discussion of renormalization in the hybrid case is more 
subtle than in the ^-derivable case, because two differ- 
ent levels of approximation for the 2PI effective potential 
are intertwined, each of which comes with its own set of 
countertcrms. First, the Hartree-Fock effective poten- 
tial, from which the hybrid gap equations are deduced, 
is obtained from Eq. (|36[) after removing the setting-sun 
sum integrals, making the replacements m-2 mo and 



A 



(A,B) 



-> A 



(A,B) 



and taking A4 as given by Eq. ([32 

(A,B) 



The parameters m-2 and X\ ' are taken equal to mo 
and Aq j4 ' B '* because, in the Hartree approximation and at 
= 0, M 2 = M 2 and V = V. To understand why fl32J} is 
the relevant choice for A4 , one has to recompute V in the 
Hartree approximation, along the lines of Appendix [X] It 
follows in particular that the gap equations in the hybrid 
approximation read 



M 2 



. (A+2B) 
,2 1 A \a2 



'0 



A 



6N 

((N-1)A) 




6N 



-T[G T ] 



(38) 



and 



A 



A o 
6N 

((iV-l)A+2B) 



T[G t 



6N 



-T[G T ] , 



(39) 



which are obtained equivalently from Eqs. ([BJ and (J7|) by 
disregarding the momentum dependent pieces and mak- 
ing the replacements X^' B ^ — > Xq A ' B \ Second, since the 
two-loop 2PI effective potential is evaluated for a dif- 
ferent propagator than in the two-loop case, the bare 
parameters needed to renormalize the hybrid effective 
potential need not be the same as those derived in the 
previous section. In fact, it can be immediately seen 

that mo, m,2, Xq A ' B \ X^' 3 ^ are the same as they are also 
needed to renormalize the gap and curvature masses at 
4> = which remain unchanged. However, as mentioned 
above, the field equation receives additional contributions 
and therefore A4 is changed. Another point of view is that 
the four-point function V^, = o is modified as compared to 
the two-loop ^-derivable case. After some calculation 



7 There was a factor of 1/2 missing in Eq. (22) of Ref. [| and as 
a consequence there should be a factor of 2 in front of the two 
terms of the last line of Eq. (31) of that reference. 



whose details are gathered in Appendix [B] and upon im- 
posing the renormalization condition V^=o,t, = A*, we 
arrive at 



A 4 = -2A. + 2A^ + ^-A - 4+2 ^ (40) 

which gives A4 = — 2A* + 6A21 — 3Ao when TV = 1. 

One of the nice features of the hybrid approximation is 
that, since the self-energies are momentum independent, 
the divergent part of the various sum integrals involved 
in the calculation can be determined analytically. 
In this way, one can check explicitly that the above 
countcrterms renormalize the gap and field equations 
as well as the effective potential and explicitly finite 
expressions can be obtained for them. We will now show 
in detail how to derive the finite gap equations and 
then sketch the derivation of the finite hybrid effective 
potential which in turn leads to a finite field equation 
by differentiation. 



B. Explicit renormalization 

Recall first how the renormalization of the gap equa- 
tion for <fr = works. We have seen in this case that 
everything boils down to a single equation, Eq. (|13[) , sim- 
ilar to the gap equation in the case N = 1. After using 
the value of m§ , one obtains 

[T[G^]-%[G+]}. (41) 



A 



(NA+2B) 



1 

6N 



By using the techniques developed in Ref. [|| or by per- 
forming an explicit calculation, it is easily seen that 
the remaining divergence in the right-hand side is noth- 
ing but -(M| =0 - m2)£*[G*](0)/6AT. Subtracting this 
contribution from both sides of the equation and using 
Eq. (|26j) . we end up with 

, N + 2 



M| =0 



67V 



-A*7f[G0 = 



(42) 



where we have introduced the finite combination 



T F [G] = T[G] - %[G ! *] + (A/ 2 - ml)B4G*}(0) . (43) 

In order to generalize these manipulations to the case 
4> 7^ 0, we note that what matters when <j> = is that the 
combination of masses appearing in the left-hand side 
of Eq. (|4Tj) is exactly the same as the combination of 
tadpoles in the right-hand side and also that the combi- 
nations of bare couplings Xq A ^ and Xq B ^ are precisely the 
ones given in Eq. (|26[) whose inverse is finite up to a bub- 
ble diagram with the appropriate prcfactor. If we were 
able to find linear combinations of the masses Ml and 
Mf. involving the same linear combinations of the cor- 
responding tadpoles, we could apply the previous proce- 
dure twice. Now, if we write the system of gap equations 
as 



M 2 



( r[G L ] 

I T[G T ] 



(44) 
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we see that if P denotes the matrix that diagonalizes the 
system, we obtain 



Ml 



A 
n 



P ( T[G h ] 
{ T[G T ] 



P 



(45) 



and thus P provides the sought-after combinations. 
These are found to be M[ + (N ~ 1)M| and M[ - M| 
and we note that the corresponding equations not only 
involve by construction the same combinations of tad- 
poles, that is T[G L ] + (N— 1)T[G T ] and T[G L ] - T[G T ], 
but also that they involve respectively the combinations 
Xq NA+2B ^ and Aq 25 '' which are those whose inverse is fi- 
nite up to a bubble diagram with the appropriate pref- 
actors, see Eqs. (|2"5"j) and (f2"u) . Wc can now apply twice 
the procedure used for the case <j> = and, after switch- 
ing back to longitudinal and transverse components, we 
finally end up with the equations 



Ml 



and 



Ml = ml 



Tf[G l 



2N 1 



N - 1 
6JV 



N + 1 
6N 



A.Tf[G t ], (46) 



A,T F [G T ], (47) 



which are both finite. A similar approach has been 
used in Ref. in the case of a theory with two scalar 
fields, not related to each other by 0(2) symmetry. 
Surprisingly, the author did not use this approach in the 
case of the O(N) model in Ref. [2(j. This is probably 
related to the fact that he was not considering multiple 
bare couplings as we do here, see also the discussion 



below. 

To sketch the renormalization of the two-loop hybrid 
effective potential, let us consider the case N = 1 first. 
The trick is to express the hybrid effective potential in 
terms of the Hartree-Fock effective potential 

7hW = 7o(mi) + vAtf + + h [r[ G] + 2 2 ]r[G] 



\j [inG- 1 - In G^ + (Q 2 



+ m 2 )G-i\ 



(48) 



which we know how to renormalize, see Ref. [15j |. Wc 
have 



l{4>) = 7h(0) + 



m 2 - m ,2 , ^21 - A 
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(49) 



where we have used the fact that in the Hartree-Fock 
approximation m 2 is equal to m , A 2 is equal to Ao and 
A4 is given by Eq. (p?2"j) instead of Eq. (|4"0"|). so that the 
difference accounts for the <fi 4 term above. Using the ex- 
pressions for m 2 and m , together with the gap equation 
at TV = 1, this term cancels and we arrive at 



7(0) = 7h(0) + y^ C [ (5 ' G *]' 



(50) 



with 



c[S - G - 1 ^(^- 1 ) (fi2 

= 7f[G]B*[G*](0)- 



ml) + 



^A 2n l 



[T[G]-%[G i ]]--[S[G]-S it [G i ]] 



1 



S[G]-S+[G*]-{M 2 -mi) 



dm? 



(51) 



The second line has been obtained by using the explicit 
expressions for A21 and <5A 2n i and shows that the determi- 
nation of C [G, G+] relies essentially on the determination 
of S[G]. An explicit proof of the finiteness of C[G, G*] is 
given in Appendix [Bj This concludes the proof that the 
hybrid potential is finite in the case N = 1. We men- 
tion that a finite expression of C[G, G*], which can be 
used for the numerical evaluation of the effective poten- 
tial, was obtained within dimensional regularization in 
Ref. 0, see Eq. (Bll) there. 

Similar considerations for arbitrary N lead to 



7(0) = 7h(0) + 



367V 



Cn[Gl, Gt, G*] 



(52) 



where 



C N [G l , G T , G*] = (N + 8)C [G L , G*] 



+(IV-1)C[G L ,G T ,G*], (53) 
with C given in Eq. (|51|) and 

1 



C[G L ,G T ,G,] = 2T F [G T ]^[G,](0) - 



35 [Gl; Gt; Gt] 



-5[G l ] - 25*[G*] - 2(M| 



m2) ^[G,] 



(54) 



As it was the case for C, it is possible to show that C is 
finite and we refer to Appendix [Bl for the details. 

Thus, it remains to be shown that the Hartree-Fock 
potential is renormalized for arbitrary N. In fact we 
expect it to be finite up to a temperature and field inde- 
pendent divergent constant. For this reason, we consider 
instead the subtracted potential Aj(cj)) = — 7*(0). 
In the case N = 1, one possibility is to rewrite the ef- 
fective potential in terms of the combination <f> 2 + T[G]. 
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We complete a square of the form oc (<fi 2 + T[G]) 2 and 
gather the terms proportional to the bare mass into 
m o( ( t> 2 + 7"[Cx|)i then we use the gap equation to write 
4> 2 + T[G] = 2(M 2 -m 2 ] )/X and the expression for 1/A , 
which can be read off from Eq. (|26|) . Performing these 
steps we end up with the subtracted effective potential 

A7h(^) = + l(U[G] - M 2 T F [G]) 



M 4 - mi 
2K 



(55) 



where we have introduced the subtracted logarithmic 
sum integral 

£ F [G] ee 2[ 7o (m*) - To(m*)] + f [^G' 1 - lnG" 1 ] 

Jq 

- (M 2 - ml)%[GA + \(M 2 - ml) 2 B4G*}(0) , 



(56) 

which can be checked to be finite. It remains to be shown 
that the combination A4 — 3Ao is finite. From Eq. (|32p . 
we have A4 — 3Ao = — 2A*, which concludes the proof in 
the one-component case. 

The extension to TV ^ 1 is rendered difficult by the 
presence of terms of the form T[Gl]T[Gt] which cou- 
ple longitudinal and transverse components. However, 
if one expresses the Hartree-Fock potential in terms of 
the diagonalizing combinations obtained above, namely 
Gl + (N — 1)Gt and Gl — Gt, one checks that such types 
of coupled terms disappear. Moreover the combinations 
of bare couplings which come with such decoupled combi- 
nations arc again precisely those for which we have simple 
expressions given by Eqs. (f2~5j) and (|2l)|) . We can thus re- 
peat twice the standard procedure for the case N = 1. 
After switching back to longitudinal and transverse com- 
ponents, we finally end up with the rcnormalized cxpres- 



A 7 h(0) 



A*0' 



4 



12N 

~(N + 2)\ ir 



-(£ F [G L ]-A^T F [G L ]) 



N+l 



(Mi mi) + 



N -I 



(£f[G t ] 



3{N- 1) 



M 2 T F [G T }) 
N - 1 



(M 2 M 2 - ml) 



(57) 



In the hybrid approximation we shall not use the field 
equation, but search for the minimum of the effective 
potential ([57]). as explained in Sec. IIV1 nevertheless, 
for completeness, we give its renormalized form in Ap- 
pendix [5] 



C. Comparison to other approaches 

To close this section, let us compare our renormal- 
ization procedure to other approaches followed in the 
literature. Since most of these approaches concern the 
Hartree-Fock approximation, we focus on the latter for 
which we have given the renormalized effective potential 
in Eq. (|57p and the renormalized gap equations in (|46[1 - 
(|47|) . The finite field equation can be obtained by plug- 
ging Eq. (|38p into the Hartree-Fock bare field equation 
to yield 



6N 



= Mr 



where we have also used Eq. (|32|) . 



J> 2 = M 2 -^ 2 , (58) 



3N 



The renormalization of the Hartree-Fock approxima- 
tion was investigated for instance in Ref. 0] where two 
different regularization schemes, cutoff and dimensional 
rcgularization, were used together with the correspond- 
ing "renormalization" schemes, named respectively "cut- 
off scheme" (CO) and "counterterm scheme" (CT) and 



leading surprisingly to different results. In fact the CO 
scheme is not really a renormalization scheme since the 
authors explain that there is no way to send A to infinity 
and the equations need to be considered at finite A, A be- 
ing an additional parameter of the model. The drawback 
of such an approach is that certain obstructions appear 
in parameter space, in particular in the chiral limit. In 
contrast, the CT scheme removes the divergences and 
the continuum limit can be considered, with no obstruc- 
tion in the chiral limit. This seems contradictory since 
one could expect that physical results should not depend 
on the regularization method used. Moreover, the CT 
scheme was not given a real justification in Ref. Q and 
it was not clear how to generalize it to higher order trun- 
cations. The renormalization that we use in this work 
clarifies these issues. As we now explain, it gives a justi- 
fication to the CT scheme of Ref. Q , it is generalizable to 
an arbitrary level of truncation and it allows us to modify 
the CO scheme in such a way that it becomes formally 
identical to the CT scheme. In particular it presents no 
obstructions in the chiral limit. 



If we have a closer look at Eqs. (TIB)) and (jT7|) for in- 
stance, we notice that, except for the fact that the sub- 
tractions are made at a finite temperature T*, our renor- 
malized gap equations have structurally the same form as 
those of the CT scheme of Ref. Q . We thus see that one 
way to justify this scheme is to admit the need for mul- 



10 



tiple bare parameter^ which need to be fixed by appro- 
priate renormalization conditions, supplemented by con- 
sistency conditions. Unlike what is stated in Ref. |8(, our 
interpretation shows that the CT scheme does not involve 
temperature dependent counterterms since the countert- 
erms depend only on the renormalization scale T* but not 
on the self-consistent mass M 2 . Moreover, these consid- 
erations are sufficiently general to be extendable to higher 
order approximations or to apply to any regularization, 
with similar results in the continuum limit. In particular, 
we can define a CO scheme for which the renormalized 
equations are PB"]) and (|77|) with integrals cut off at some 
scale A. In this scheme the cutoff A can be sent to in- 
finity (as mentioned above, this is a peculiarity of lower 
order approximations) and no obstructions appear in the 
chiral limit. In fact the problems with the CO scheme in 
Ref. @ can all be identified with the use of one single bare 
coupling, instead of multiple ones as we propose here. To 
illustrate this, let us revisit one of the obstructions raised 
in Ref. [1] and see how it is lifted within our approach. 
In the CO scheme of Ref. [|[ , the gap equations at finite 
A are written using a single bare coupling. This amounts 
to replacing A ^ and by Aq in the bare gap equa- 



tions ((38)) and ([39)) . Similarly the bare field equation is 
written with the same coupling Ao everywhere and reads 
in the Hartree-Fock approximation: 



= m 



+ 2N T[Gh] 



(iV-l)Ao 
6N 



Ml 



(59) 



1 



Ap_ 

3N' 

to be compared to Eq. ([551) . Writing the difference of the 
two gap equations at T = 0, setting the pion mass to zero 
and using the field equation (|59")) . one arrives then at 

[■T=0 1 r T=0 

° = / a* 

JQ<A W 

whose solutions are cither m a = or A = 0, both absurd. 
This is the conclusion reached in Ref. Q. In contrast, 
within our scheme, if we subtract the renormalized gap 
equations (|46)) - ([47l) and use renormalized field equation 
([55| . we obtain 

j-T=0 -] r T=0 

<A Q 2 + ml J Q<A Q 2 



Q<A V 



T, 



1 



(61) 



,<a (Q 2 + ™ 2 ) 2 ' 

which admits a nonzero solutior@ for m a , pretty in- 
sensitive to the large values of A because Eq. (|6"Tj) is 



8 As it is explained in Ref. 1151 , the need for multiply defined bare 
parameters is a truncation artifact and, the consistency condi- 
tions are such that, if one increases the order of the truncation 
the differences between the various bare parameters, should be- 
come smaller and smaller, at least formally. 

9 As a function of m^, the right-hand side of the equation starts 
at when = 0, and decreases first before growing linearly as 
rr?„ — > co. 



renormalized. 

Our approach differs also from that used by Amelino- 
Camelia and Pi in Refs. [13, HH where only one bare 
coupling was used. If we were to use only one bare 
coupling, the first term of Eq. ([55)) would be — (Ao/12)0 4 
in place of — (A*/12)0 4 . According to Amelino-Camelia 
this term does not spoil the renormalizability because 
Ao, albeit being a bare parameter, approaches 0~ as 
A — > oo. However, as already discussed above, the 
possibility to send the cutoff to infinity is a peculiarity 
of the lowest order approximations, not shared by higher 
order ones where we expect physical quantities not to 
be defined above the Landau scale. It is thus more 
satisfactory to implement a renormalization scheme in 
which the results are already pretty much insensitive 
to the cutoff below the Landau scale. This is achieved 
by our scheme if the Landau scale is not too close to 
the physical scales because our results show a "plateau" 
behavior below the Landau scale, whereas in the scheme 
by Amelino-Camelia there remains a logarithmic sensi- 
tivity from the term — (Ao/12)0 4 . One could argue that 
the existence of a plateau is related to the existence of 
a continuum limit, a notion that does not make sense 
at higher orders of approximation. Still, as already 
mentioned above, if the parameters are such that the 
Landau scale is much larger than the relevant physical 
scales, there is an intermediate regime where this notion 
can be considered in a somewhat generalized accepta- 
tion: quantities renormalized according to our scheme 
will still show a plateau behavior for values of the cutoff 
below the Landau scale. These considerations can be 
made more quantitative by using specific examples and 
will be presented elsewhere [22J . In our present two-loop 
approximation we will study the cutoff dependence 
of some physical quantities for different values of the 
parameters, that is different values of the Landau pole 
(see Figure [5]). 

Let us finally mention that certain works disregard 
renormalization by arguing that one is only interested in 
thermal effects and thus that "vacuum" fluctuations can 
be neglected, see for instance Ref. Q. It is worth men- 
tioning however that, in a self-consistent context such as 
the 2PI formalism, the masses or self-energies that enter 
these vacuum fluctuations depend on the temperature. 
Neglecting them is then not completely justified and can 
lead to neglecting an important piece of the thermal con- 
tribution. This can be tested by using the exact limits 
of certain models/theories such as the limit of a large 
number of flavors in QED/QCD, see Ref. [23j . 



IV. NUMERICAL METHOD 

Before discussing our results in the next section, let us 
give a brief overview of the numerical methods that we 
used to solve the equations and compute various quanti- 
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ties of interest. 

In the two-loop case we take advantage of the fact 
that all momentum dependent sum integrals are convolu- 
tions and compute them by means of discrete fast Fourier 
transform algorithms (DST and DCT as described in 
Refs. Sil) 

using a 3D cutoff A. We exploit the rotation 
symmetry of the propagators to reduce our discretization 
to a two-dimensional N T x N s lattice containing 7V r — 1 
positive Matsubara frequencies in addition to the static 
mode L) n = and N s moduli of the 3D momentum, the 
smallest available being the lattice spacing in momentum 
space Ak = A/N s . Moreover, since the leading asymp- 
totic behavior of G{Q) is exactly 1/Q 2 in the approxi- 
mation at hand, we can increase the rate of convergence 
of the Matsubara sums and of the convolutions by sub- 
tracting first the leading (free-type) asymptotic behavior 
of the various summands/integrands. These subtracted 
sum integrals involve free- type propagators, as it is also 
the case for all the sum integrals encountered in the hy- 
brid approximation, and therefore can be computed al- 
most exactly. In practice this means that the Matsubara 
sum is performed exactly and the momentum integral is 



computed numerically using accurate adaptive integra- 
tion routines of the GNU Scientific Library (GSL) j25| . 
For more details on the numerical aspects, we refer to 
our previous work Q and adopt the notations used in 
its Sec. V. In the remainder of this section, we describe 
some of the most important aspects, in particular the 
new features that appear in the case N ^ 1. 

A. Increasing the rate of convergence of the sum 
integrals 

After using the expressions for ttiq, m\ and 5\^[ B \ 
which can be read off from Eqs. (pi) . ([22]). ([27]) and ([2"9]t. 
it is straightforward to apply the procedure described in 
Sec. V.B of Ref. @ to render the longitudinal gap equa- 
tion (|6|), the gap equation at <f> = (fT3|) . and the expres- 
sion of the curvature at <fi = (fl5|) in a form suitable for 
numerical computations, because they contain the same 
types of sum integrals as those in Ref. [ij. This is true 
also for the subtracted effective potential, defined using 
Eq. (|36j) as A"f((f>) = 7(0) — 7*(0), when it is written using 
the expression J7 /5<j) appearing in the field equation^ as 



A7(^)=Ar( 7o (m*,A)- 7 S(m* ! A))+ £ f / \^Gr\Q) -\aG~\Q)- {ffi(Q) -m$Gi(Q)] 
A 4 



i=T,L 



+ 



-2% 
12AT 



24N 



4 _ H 

2 2AN 



T[SG L ] + (N- 1)T[<SG t ] + N6T[G*] 



(T[SG L ] + 5T[G4) 2 + (N- 1)(T[6G T ] + «JT[G*]) 5 



(62) 



where cl = 1, ct = N — 1, 7o(m + ,A) is the inte- 
gral in Eq. (|37p calculated with a cutoff A and at a 
temperature and we used the shorthand notations 
SG h/T = G L/T - G* and ST[G*} = T[G*] - %[G*\. The 
integrals are evaluated as shown in Eq. (131) of Ref. Q. 
However, for the transverse gap equation and the field 
equation itself (|14[) . we need to compute two sum inte- 
grals which were not encountered in our previous work. 
The first is the bubble sum integral of Eq. ((TJ) , which is 
rewritten as 

B[G L ; G T ](K) = B[G*}(K) + [ G L (Q)5G T (K - Q) 

JQ 

+ f SG h (Q)G i ,(K - Q), (63) 

JQ 



In the presence of an external field h, S"f/S(j> is the expression 
appearing on the right-hand side of Eq. Q , but multiplied by 
<j> and with h</> subtracted from it. 



where 6Gl/t decrease faster in the UV than Gl/t, hence 
reducing the error of the corresponding sum integrals, as 
compared to that of the sum integral £>[Gl;Gt], while 
the first term involves only the free-type propagator G* 
and can be computed almost exactly. The discretized 
form of Eq. ([63f used in the numerics can be easily given 
in terms of the discrete version of the convolution defined 
in Eq. (114) of Ref. @. The second new sum integral is 
decomposed as 

5[G l ;Gt;G t ]= / G L (Q)(B[G T ](Q) - B[G*](Q)) 

JQ 

+ [ 5G L (Q)B[G,](Q) + S[G^}.(64) 

JQ 

The third term and the bubble B[G*](Q) in the second 
term can be computed almost exactly. The summand in 
the second term decreases faster than the original one, 
Gl(Q)B[Gt](Q) , as it is the case with the difference of 
bubbles in the first term, which can be rewritten as a con- 
volution using Eq. (121) of Ref. @. Then, the discretized 
form of 5 [Gl ; Gt; Gt] can be readily written using the 
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discrete version of the convolution and of the local sum 
integral defined in Eqs. (114) and (115) of Ref. @. 

B. On the solution of the equations 

In the two-loop approximation the solution of the gap 
equations ([5]) and ([7]) either at fixed </> or together with 
the field equation (fl"4"]l is obtained iteratively. In both 
cases the coupling counterterms X^' B \ \^ B%> and A4 are 
evaluated first using accelerated Matsubara sums, as ex- 
plained in Appendix C of Ref. . Then, the T-dependent 
integrals which do not depend on the solution of the equa- 
tions are evaluated using adaptive numerical integration 
routines. The quantities determined up to this point are 
unchanged during the iterative process. The process used 
to solve the coupled equations ©, ([7]) and ([T4"|) at h ^ 
is similar to that used in Ref. Q. At a given T both 
propagators are initialized with G*. The iteration starts 
with the evaluation, using the most recent G l /t> of the 
local-type sum integrals in the field equation, which is 
easily solved for it is cubic in <j). Using the obtained value 
of 4>, the propagators are updated sequentially, starting 
with Gl- First, the self-energy M*(iu) n ,k) is evaluated 
by computing the required sum integrals with the most 
recent propagators (due to the sequential update of the 
propagators there is no need to recalculate all the local- 
type sum integrals). Then, the updated propagator is 

G L (iu n , k) = [uj 2 n + k 2 + aMt + (1 - a)Mj; i0ld ]-\(65) 

where "old" refers to the propagator of the previous 
iteration, which has to be stored. The updated Gl is 
then used to update Gt in an analogous way, using the 
same a G (0, 1] parameter, which controls the speed of 
convergence of the iterative process. For large A* one 
needs a < 1 for the iteration procedure to converge at 
all, however, for small couplings the fastest convergence 
is achieved with a = 1. Besides (f>, the value of the 
propagators at the lowest available frequency and mo- 
mentum is also monitored. The iteration stops when the 
relative change of all these quantities from one iteration 
to the next is smaller than the desired accuracy (usually 
a relative change smaller than 10 was required). 

In the hybrid approximation the gap equations (|46[) 
and (|47[) are momentum independent and, therefore, 
much easier to solve compared to the full two-loop case. 
However, the field equation is complicated due to the fact 
that the propagators do not fulfill the stationarity condi- 
tions. For this reason, we evaluate instead the effective 
potential (|52[) and search for its minimum. During this 
process the vacuum parts of the sum integrals can be 
calculated analytically, while the explicitly temperature 
dependent parts can be computed almost exactly using 
adaptive numerical integration. Note that one can avoid 
the determination of Ml/t as a solution of two coupled 
equations. In the next section, we will see in Eq. (|71[) 
that it is possible to explicitly express Ml in terms of 



Mr- Plugging this expression into Eq. (|47|) yields a one- 
dimensional equation for Mt, to be solved for any <f>. 
Then, Mr, Ml at the minimum of the potential are eas- 
ily obtained with a numerical minimum finder routine, 
which chooses values of 4> and checks the value of the 
effective potential (f5"2"f evaluated with the solution Mx 
of the one-dimensional gap equation and Ml determined 
from Eq. ([TTj) . 

C. Determination of the (pseudo-) critical 
temperature and zero temperature quantities 

In the chiral limit the critical temperature T c is the 
value at which the curvature of the potential vanishes 
at (f) = 0. Since the latter is the same in both approx- 
imations considered in this work, the corresponding T c 
is also the same. Moreover since the gap equation (TT5)) 
yields a momentum independent solution, the curvature 
at vanishing field, and therefore T c , can be evaluated al- 
most exactly, using adaptive integration routines. Using 
Eqs. d^nj) and d3UJ) in Eq. ([15]) one can even obtain 

an explicitly finite equation for the curvature at vanish- 
ing field in terms of the the gap mass given by Eq. (|42|) 
and C[G0 =o , G*], defined in Eq. (f5"Tj) : 

M| =0 = M| =0 + ^±l\lC[G^GA. (66) 

The critical temperature T c is then obtained from the 
previous expression through the relation M^ = o.t c = 0. It 
is also convenient to define a temperature T c from the 
vanishing of the gap mass at (f> = 0: M^ =0 y c = 0. The 
temperature T c is the same in both approximations and 
can be given analytically, since by setting M^ = o to zero 
in Eq. (|4"2"f and introducing 

C * = m * + ^r^7~F,T=o[G ], (67) 

with Gn(Q) = 1/Q 2 , as in the one-component case in 
Ref. [ll, one obtains f c = [-72NC*/{{N + 2)A*)] 1 / 2 , if 
the parameters are such that G* < 0, otherwise it is not 
defined. We note that because the gap and curvature 
masses admit a continuum limit, so do the critical 
temperatures. Of course these continuum values are 
not directly connected with the critical temperatures of 
the systems the model could describe at low energies, 
because these are nonuniversal quantities which depend 
on the microscopic details of the particular system under 
study. To obtain them, one should rather envisage a 
first principle calculation or include sufficiently enough 
nonrenormalizablc operators in the model. 

In the physical case, i.e. at nonzero h, the pseudocrit- 
ical temperature T pc is defined through the inflection 
point of the 4>{T) curve, which is determined in both 
approximations with the same algorithm. This takes into 
account that since in the two-loop case the computation 
is time demanding, it becomes worthwhile to determine 
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the inflection point by running the code at the least 
possible number of temperature values without giving 
up the accuracy requirements. As a first step of the 
algorithm we compute 4> at five equidistant temperature 
values between T c and min(3</)(T c ), 5T*) (this proved 
always larger than T c ), where T c is the critical tempera- 
ture corresponding to the actual value of the parameters 
to 2 /T 2 and A*, but h = 0. Then, from this set of points 
we compute numerically the first and second derivatives, 
using the highest possible order of finite difference for- 
mulas for central or one-sided approximations prij which 
can be reached at a certain value of the temperature, 
given the finite number of points we have. Using the 
information that d(j)/dT has a minimum at T = T pc and 
d 2 <t>/dT 2 changes sign as it goes through T = T pc , where 
it vanishes, we can determine from our five points the 
two values of temperatures T < and T> which enclose 
the inflection point (T< < T pc < T>). Next, a rough 
estimate for the pseudocritical temperature, T ost , is 
obtained from T< and T> through a linear interpolation. 
Finally, we compute <j) at three more temperatures: 
(Test + r<)/2, (r^ + T>)/2 and by fitting the 

function f(T) = a + 6arctanc(7r — dT) to the values of 
<f> available at these temperature values and at T< and 
T>, we obtain our best estimate for the abscissa of the 
inflection point: T pc = n/d. 

The determination of the T = quantities required for 
the parametrization of the model or for the computation 
of the pressure, see below, is different in the two approx- 
imations that we consider. In the hybrid case the vac- 
uum parts of the integrals can be evaluated analytically, 
rendering the T = value of M L / T and <f> or even the 
effective potential easily accessible. On the contrary, in 
the two-loop approximation it is impossible to explicitly 
reach T — due to the use of a finite number of Matsub- 
ara frequencies, since the number of needed frequencies is 
inversely proportional to the temperature. However, this 
shortcoming of the numerical method is overcome with an 
extrapolation procedure which uses the low temperature 
data obtained by increasing the value of N T to an appro- 
priate value (see the caption of Figure [7] for an explicit 
example). The effective potential at T = is obtained 
by fitting to the low-T values a functional form based 
on the temperature dependence of the ideal gas pressure, 
g{T) = a- &T 5 / 2 cxp(-c/T). To obtain M L/T and <j> at 
T = wc use a fitting function j(T) = a — 6exp(— c/ y/T), 
which has a purely empirical motivation. 



Landau scale is large enough!"! We note that, for N = 1, 
Eq. (35) goes over into Eq. (48) of Ref. @|3 which means 
that the scale of the Landau pole obtained in the N = 1 
case for some value of the coupling is obtained in the N = 
4 case at twice that value. The value of the Landau pole 
corresponding to a given A* can be accurately estimated 
for m* <C A p using the formula 

1 -8tt 2 £W =oo [G*](0) , 



A" 



■ cxp 



(N + 2)A* 



(68) 

obtained by replacing in Eq. (|3"5j) the "thermal" part of 
the bubble integral [G*](0) with s£l =tx) [G*] (0) (we 
use the notations of Ref. @ ) . 

We also use the T c = curve, whose equation 
AJ(to*/T*) can be si mply obtained using Eq. (|67p from 
the relation C* = 0, airi 
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mlT F m =0 [Go]. 



(69) 



N + 2 

This curve can be seen in Figure [T] For points which 
are above (below) it < (C* > 0). Wc shall need 
similarly the curve T c = whose equation A^(m*/T*) is 
obtained implicitly from A/0 = o,t c =o = 0, with the renor- 
malized curvature mass given in Eq. (|66[) . that is 

"A*Ct c =o[G^=o,t c =0) G-k]- (70) 



— Mj,_ T o 
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We note that the solution of the gap equation at vanish- 
ing temperature (and for N — I even at arbitrary value 
of the field) can be obtained in closed form in terms of 
the two real branches of the Lambert function W. In 
Appendix [B] we provide the solution at vanishing field 
and temperature, which can be used in (|70[) to obtain 
A^(to*/T + ) numerically. 

The last characteristic curve is needed in the hybrid 
case. In the Hartree-Fock approximation applied to the 
one-component case it was already observed in Ref. (l5j 
that there is a temperature dependent critical value of 
the field, <j> c {T), satisfying (j> c {T) < 4>(T), such that for 
smaller values of the field the gap equation does not ad- 
mit a physical solution. Wc investigate now the exis- 
tence of such a curve in the hybrid approximation and its 
location with respect to <j)(T). Subtracting three times 
Eq. (|47l) from Eq. (l46l) one can express Ml in terms of 
M T as 

2mi-^^KT F [G T ]. (71) 



3M£ 
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D. Characteristic curves 

In preparation for the discussion of the results in 
Sec. [V] it is convenient to define certain characteristic 
curves in the parameter space (m 2 /T 2 , A*). 

A first class of curves that we use is made of the iso- 
A p curves which allow us to determine a region where the 



11 We have treated the hybrid approximation using dimensional 
rcgularization and taking the continuum limit after proper renor- 
malization of the equations. We could have proceeded equiva- 
lently using a 3D cutoff as in the two-loop case. The iso-A p 
curves need to be understood in this context. 

12 Note that there is a factor of 1/2 missing in front of the integral. 

13 A simpler expression, AJ(m*/T*) « (72Nml /T?)/[(N + 2)(1 - 
3m*/(27rT-t))], is obtained using high-temperature expansion, 
which is reliable for m* /T* < 1 and sufficient for our purposes. 
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Expressing 7f [Gt] from the relation above and plugging 
it in Eq. (|4T|) . one obtains 

(1 - N)M% = 2to 2 ~{N + l)Ml 

+^ F (7V + 2)(0 2 +r F [G L ]). (72) 

We define 4> C {T) as the value of the field for which 
M T (T) = 0. Then, from Eq. (J72J) one obtains 

VcV y (iV + 2)A* [ L ' cV ' N + l\ 1 ' J ' 

(73) 

where M£ C (T) = A*(N + 2)(T 2 - T 2 )/(36iV) is obtained 
from Eq. (Hi}- Note that, by definition of T c , M T (T C ) 
vanishes at = 0. It follows that (f> c (T c ) = 0. The exis- 
tence of the 4> C (T) line depends, of course, on the values 
of the parameters. It is an interesting question what 
happens at zero temperature because there could exist 
a C: o curve in the (m 2 /T 2 ,A*) parameter plane along 
which <pc(0) = 4>(0) 7^ 0, and which delimits a parameter 
region where the model cannot be solved at T = in 
the hybrid approximation. This is actually the case to 
the left of the </> C) o line in the parameter region shown in 
Figure [TJ Wc have not seen any trace of a </) Ci o curve in 
the two-loop approximation. 

We mention that we could have alternatively defined 
4> C {T) as the value at which Ml(T) vanishes. In this case 
from Eq. ([721) we have M| iC (T) = -(N + 2)A*/(3(iV - 

1) )[4>1(T) + (T 2 - T 2 )/12],' which is only positive if for 
< T < T c one has 2 < (T 2 - T 2 )/12. This means that 
one needs C* < 0. Then similarly as in Appendix IB1 (one 
only has to change b± given in Eq. (|B18|) to 6* — 1/2) the 
T = solution of Eq. (j7ip at Ml = given in terms of the 
Lambert function is bigger than 2A° st exp[247r 2 7V/((7V + 

2) A*) — 1]. The size of this solution matches the size of 
the Mt, c (0) expressed from Eq. (|72)) only for very large 
A*, when the scale of the Landau pole is small. This is 
outside the region of the parameter space wc would like to 
investigate. It is easy to see, using that 7f[Gt] increases 
with T at fixed M T , that for < T < f c Eq. ([71]) admits 
only a large scale solution, because the right-hand side of 
Eq. ((TTj) at Mi = is negative in this temperature range 
and vanishes only at T = T c . In conclusion, in the region 
allowed by the definition of 4> c (T) based on the vanishing 
of Mr, it turned out that Ml is always positive when 
Mt is nonvanishing. As a last remark we note that in 
the hybrid approximation for the one-component case, 
where the only possible definition for tp c (T) is the one in 
terms of Ml, we can prove that the curve 4>c,o does not 
exist, because, as discussed in Appendix [Bl </> c (0) = (f>(0) 
cannot happen for (f)(0) > 0, that is for parameters for 
which the model is in its broken phase at T = 0. 

V. RESULTS 

In this section we present our numerical results on the 
phase transition and the thermodynamic properties of 



the model. As it was the case in the one-component 
scalar model studied in Ref. @, we shall find that in 
the chiral limit, the transition, when it occurs, is of the 
second order type. We shall show this explicitly by mon- 
itoring the variation of the order parameter. We shall 
also determine some critical exponents as well as thermo- 
dynamical obscrvables. Before doing so, we discuss the 
physical paramctrization of the model, relevant, when 
— 4, for the discussion of light meson properties. 

A. The parametrization of the model 

The renormalized O(N) model has three parameters 
to 2 , A* and h (h = in the chiral limit) and a renormal- 
ization scale T*. Being the solution of the gap equations 
at (j) = and T = T*, m 2 is positive, and since we want 
the bare couplings to be positive, we need to restrict to 
A* > (in addition to A < A p ). Not all the 4-uplcs 
(to*, A*, h, T+) correspond to different physical systems. 
First of all, renormalization group invariancc implies that 
given two values for the renormalization scale T*, there 
exists a renormalization group transformation that maps 
two sets of values for m* and A* in such a way that the 
physical predictions are the same. This is rigorously true 
in the exact theory where no approximation is consid- 
ered but it needs not be the case in a given truncation of 
the ^-derivable potential and the T* dependence of the 
physical results needs to be investigated, see below. 

Another source of redundancy is provided by dimen- 
sional analysis, since knowing the values of the physical 
observables of a system represented by (to*, A*, ft., T*), 
one can very easily deduce the values of the same phys- 
ical observables for a rescaled system represented by 
(am t ,A t ,a 3 /i,aT t ), where all dimensionful quantities 
are rescaled by a to the appropriate power. In contrast 
to renormalization group invariancc, this redundancy is 
present at any level of truncation and it is therefore con- 
venient to get rid of it by working exclusively with di- 
mensionless parameters. For instance, below we will be 
interested in the value of the order parameter at T = 0, 
which is a function 4>q = ^(m 2 , A*, h; T+) of mass dimen- 
sion one (the label emphasizes that the given quantity is 
computed at T = 0). Using simple dimensional analysis, 
we deduce that 

^o/T* = Mml/Tl A*, h/T*; 1) . (74) 

Similar expressions can be obtained for the rescaled 
curvature masses Ml.o/T* and Mt,o/?* that are also 
needed below. The use of rescaled variables to 2 /T 2 
and h/T^ as parameters, is more suitable for numerical 
calculations for only dimensionless numbers are used 
and according to Eq. ([74"]) we can replace by 1 in the 
numerical code. 

In principle, the parameters can be fixed by equating 
quantities computed at zero temperature with their ex- 
perimental values. Our choice is to relate 4>q with the 
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pion decay constant / w and the curvature masses Mr,o 
and Mi^o with the mass of the pion and sigma parti- 
cles, m n and m CT respectively. We decided to use those 
masses for they reflect the best symmetry of the the- 
ory whereas, as discussed in Sec. Ill A| the gap masses 
violate the Ward identities associated to the O(N) sym- 
metry, e.g. Goldstone's theorem. However, the choice 
of curvature masses for paramctrization is questionable, 
since usually the measured physical masses are the pole 
masses. In this work, we do not have access to the spec- 
tral functions and therefore we assume implicitly that the 
pole masses are not so far from the curvature modes (this 
of course would deserve further investigation) . 

One way to proceed would be to choose a value of 
and equate 4>o in Eq. ([71)1 to / w : 

4> a = T+Mm 2 jTl A., h/T*; 1) = f n (75) 

and similarly for m a and m w . This would define a point 
in the parameter space (m+/T+, A*, h/T+). By chang- 
ing the value of T+ without changing the values of 
m a and m w , we would then follow a line of constant 
physics. One difficulty with this approach is that our 
renormalization procedure requires in the chiral limit 
the temperature T+ to be necessarily in the symmet- 
ric phase and thus for a given set of physical values of 
/tt, vn a and m„ there is a minimal possible value for T* 
which we do not know a priori. Another difficulty is 
that the sig ma mass is not known exactly, as according 
to Ref. [23 m a S (400, 550) MeV and, based on large- 
TV studies (28l . [29| . one may have concerns whether in 
our approximation Ml,o turns out to be large enough 
Hence, instead of trying to fix the parameters by picking 
up some arbitrary value for the sigma mass in the range 
given above, our procedure is to scan an appropriately 
large part of the space (mJ/T*, A*, h/T+) and determine 
at each point 4>o/T+ using Eq. (fT4")l and Mt,o and Ml,o 
from similar relations. At each point of the investigated 
parameter space we require 0o = 93 MeV, which fixes T* 
according to Eq. (|75l) and allows us to determine Mt.o 
and Ml,o- We then keep only those points which satisfy 
M T = 138 ± 1.38 MeV and allow for the decay of the 
sigma particle into two pions by requiring Ml,o > 2Mt,o- 
A one percent tolerance is allowed in the value of Mt,o 
in order to guarantee a sufficient number of points, even 
when the parameter space is not densely sampled. In 
the chiral limit, there is no constraint on Mt.Oj because 
this vanishes due to Goldstone's theorem (see the dis- 
cussion below Eq. (|17|) ) , and hence the constraint on the 
sigma mass is lifted as well. Another difference is that 



A maximal value of the sigma pole mass was observed in these 
studies. This can be seen in Figure 2 of Ref. [2H , and a formula 
determining the maximal value was derived in Ref. [29j . The 
renormalization scale used to fix the coupling constant differs in 
the two references. 



the value for the pion decay constant in the chiral limit 
is f%=° = 88 MeV [U instead of f n = 93 MeV used at 
h ^ 0. 

Since by construction all points that we keep are such 
that 4>o and Mj.o are fixed, the iso-AfL,o curves are "lines 
of constant physics" . We use quotation marks because, 
as already mentioned, in a given truncation, we expect 
physical quantities to vary slightly as we move along 
such a line, that is as we change T* for fixed f v , m a and 
TO w Along such a line we can determine in particular 
T pc (T c at h = 0) from the inflection point of the 4>(T) 
curve as described in Sec. QV] and plot its dependence 
with respect to T*. We can apply the same strategy 
for any other physical quantity and study its physical 
dependence, see Sec. IV Dl for a discussion concerning the 
pressure. Note finally that, even though all points in 
parameter space correspond by construction to a given 
value of </>o , we could access other values of </>o (if our 
model would apply to other physical systems) by using 
dimensional analysis and changing the corresponding 
value of T+. Of course, all other dimensionful quantities 
would be scaled by the same quantity. 

The result of the paramctrization in the chiral limit 
is shown in Figure [1] in an almond-shaped range of the 
parameter space encountered already in Refs. The 
Ap/T* = 50 curve can be easily obtained from Eq. (|35p or 
Eq. ([S8]). the f c = curve is given by Eq. (|69|). while the 
T c = curve is obtained by solving Eq. ([70)1 using (|B19I) . 
The points investigated in the two-loop case are shown 
in Figure [T] by squares in order to distinguish them from 
those used in the hybrid approximation which populate 
more densely the studied region and appear in the form 
of vertical lines. 

In the hybrid approximation, the region to the left 
of the <p c fi curve is excluded at h = because, as 
discussed in Sec. IIVDI the model cannot be solved at 
T = in that region. Actually, the presence of this 
line, along which Mt,o = 0, invalidates the use of the 
hybrid approximation in the chiral limit in a relatively 
large region of the parameter space, the grey region 
of Figure Q] This is because as one enters this region, 
by decreasing for example ml/T 2 at fixed A*, Ml,o 
increases very abruptly. Such a huge sensitivity to 
the parameters alone raises suspicion concerning the 
applicability of the approximation, but in our case 
one can check explicitly that the results of the hybrid 
approximation deviate in this case from those obtained 
in the full two-loop approximation. The right boundary 
of this region is given by the points where the relative 
change of Ml,o compared to the two-loop approximation 
equals 3%. Apart from this excluded region, the results 
obtained in the two approximations are very close to 



15 Also, even though in the exact theory, the lines of constant 
physics should have a constant h, this does not need to be the 
case in a given truncation. 
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FIG. 1. Parametrization in the chiral limit (h — 0). The 
scanned region is bounded by the A p /T* = 50 (upper) and 
T c = (lower) curves. The <f> c ,o curve is only present in the 
hybrid approximation, in which case the grey region is ex- 
cluded for a reason explained in the text. The points which 
form vertical lines are obtained in the hybrid approximation, 
while the squares denote the solution of the full two-loop ap- 
proximation. The iso-Tc and the iso-ML,o curves are obtained 
in the hybrid case. The palette shows the value of the renor- 
malization scale T*. The inset shows the variation of T c with 
T* along iso-ML,o curves. 



each other. The value of A/l.o (sigma mass) which can 
be reached is relatively low, less than 300 MeV, and the 
critical temperature is in the range [135, 190] MeV. The 
scale at which the renormalization and consistency 
conditions arc imposed varies in a relatively large inter- 
val. Once determined, it allows us to access the value of 
the Landau pole A p in physical units and one sees that, 
in the range of the parameter space where the sigma 
mass is the largest, A p > 8.5 GeV. The inset shows the 
dependence of T c on T+ along a line of constant physics. 
Interestingly, as one goes to larger values of m 1 / along 
these lines, that is as one increases T+, the dependence 
becomes linear. 

The result of the parametrization when h ^ is shown 
in Figure [2] for the hybrid approximation. Compared to 
the chiral limit we see an increase in the value of Ml,o 
and of the pseudo critical transition temperature T pc and 
a significant decrease in the value of the renormalization 
scale T*. For fixed mJ/T*, larger values of Ml,o can be 
achieved for higher A*, that is allowing the Landau pole 
to come closer to the physical scales. We note that a 
similar figure could be obtained in the two-loop approx- 
imation, but with a significantly increased numerical 
effort. In the hybrid case the code is much faster than in 
the two-loop case and hence one can run it for a much 
larger number of points of the parameter space. We 
have tested on a good number of points of the scanned 
region, even those not satisfying Ml > 2Mt, that for 
a given set of the parameters the two-loop results for 



Ml,0i Mt,o an( i ^pc are within 3% of the values obtained 
in the hybrid approximation. This is shown in Figure [31 
where the general tendency is that at fixed m+/T+ both 
A* and h/T+ tend to increase the difference, so that the 
largest difference is obtained at the largest A* and h/T+, 
and that this largest difference decreases with increasing 
ml/Tl 

Figure 0] shows the variation of the pseudocritical tem- 
perature with the renormalization scale T* determined 
during parametrization in the physical case and in the 
hybrid approximation. The lines of the figure belong 
to different curves in the (mJ/T^, A*, h/T+) parameter 
space selected by different values of Ml.Oj each of them 
being a line of constant physics. One sees that the T* 
dependence is less than 10%. In units of T+, both T c 
and T pc decrease for increasing T* and for large values 
of T* one can fit (up to possible logs) a + b/x on T c /T* 
and Tpc/T*. In both cases b > 0, but in the chiral limit 
a > 0, while for h ^ one has a < 0, which accounts for 
the increase of T c and decrease of T pc seen in Figure Q] 
and in Figure [5] for a given line of constant physics and 
for large T+. We expect \a\ to diminish as we increase 
the order of truncation. 



B. On the sigma mass 

The parametrization reveals that there is a large re- 
gion of the parameter space where a separation of scale 
occurs in the sense that the physical scales are much lower 
than the cutoff, which in turn is much smaller than the 
scale of the Landau pole A p . In this case the solution 
of the model is practically insensitive to the cutoff used, 
as it was the case for N = 1 in Ref. H, where the cut- 
off dependence was thoroughly investigated. We have 
also seen that the value of the zero temperature sigma 
mass defined through A/l,o increases with increasing A*. 
We have reached values of sigma masses which are larger 
than the maximal value of the sigma pole mass found 
within the large- N approximation in Ref. (28j . which in 
the chiral limit is 7710- sa 328 MeV obtained for a coupling 
A « 311 and a renormalization scale of Mo Ra 334 MeV 
and m a ss 362 MeV in the h ^ case, obtained for 
A sa 386 and M sa 381 MeV. The scale of the Lan- 
dau pole in these cases is approximately 1853 MeV and 
1150 MeV, respectively. In Ref. (29|, where the renor- 
malization scale and the value of the coupling were cho- 
sen differently, a higher value of the sigma pole mass of 
around 433 MeV was reported. However, in that case, 
the scale of the Landau pole was only 720 MeV which 
prevented calculations above T sa 50 MeV. 

We investigate now what happens in our case with the 
scale of the Landau pole, which in view of Eq. ([35]) de- 
creases with A* when all the other parameters are kept 
fixed, if a more realistic parametrization of the model is 
required, in which m a £ [440, 470] MeV to conform to 
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FIG. 2. Parametrization at h 7^ in the hybrid case. The location of the investigated points relative to the characteristic 
curves is indicated with smaller size points in the (A*,m*/T^) plane of the figure in the left panel. The points having bigger 
size indicates the parameters used in Figure [3] to compare the result of the two-loop and hybrid approximations. The smaller 
size points satisfy the two criteria Mt,o = 138 ± 1.38 MeV and Ml,o > 2Mt,o- The value of the renormalization scale is 
indicated on the figure in the right panel. The two palettes show the values of Ml,o and T pc , respectively. 
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FIG. 3. Comparison between the two-loop and hybrid approx- 
imations based on the relative change of the T — curvature 
masses used for parametrization and the pseudocritical tem- 
perature (AQ[%] = 100(Qtwo-ioo P /Qhybrid - 1))- Different 
point types denote different values of m*/T* and for each 
corresponding value of A* there are points at three values of 
h/T't : 0.5, 1.5, 2.5, which generally increase from bottom to 
top. 



recent dispersive analyses of more precise irir scattering 
data (see Ref . [27| and for a recent review Ref . [3l| , in par- 
ticular its Figure 3). To this end, we have chosen differ- 
ent values of mJ/T^ and increased the value of A* in the 
range between the A p /T^ = 50 and A p /T^ = 20 curves 
of the (mJ/T^, A*)-planc, shown in Figure [T] It turns 
out that in the two-loop approximation it is possible to 



FIG. 4. The dependence of the pseudocritical temperature 
T pc on the renormalization scale T* in the hybrid approxima- 
tion at ft 7^ along different lines of constant physics specified 
by the value of Ml,o- 



reach with the parametrization procedure described in 
the previous subsection values of the Ml,o m the de- 
sired range. For instance, we obtain Ml,o ~ 465 MeV 
and T+ sa 167 MeV for ml/T? = 0.04, 'h/T? = 0.38, 
A* = 19.2, and A/ L , « 445 MeV and 171 MeV 

for mllTl = 0.124, 'h/T* = 0.355, A* = 32.476. In these 
cases the scale of the Landau pole remained at least seven 
times larger than the largest mass scale given by Ml , that 
is A p ~ 3.4 GeV. 

The interesting question is whether the scale of the 
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FIG. 5. The cutoff dependence of the relative change of <f> (upper raw) and that of Ml (lower raw) obtained in the two-loop 
approximation at different temperatures: T — 0, T = T c and T = 2T C . The different parameter sets are: (a) m+/T+ = 0.124, 
A* = 22.28, hjT'i = 1.775 for which M L ,o « 280 MeV, A p « 186 GeV and T* « 101 MeV; (b) m 2 /T 2 = 0.04, A* = 17.39, 
h/r* = 0.6, for which Ml,o « 360 MeV, A p « 16.2 GeV and T* « 146 MeV; (c) mJ/T 2 = 0.04, A* = 19.2, h/T? = 0.38 for 
which M L ,o « 465 MeV, A p « 3.35 GeV and « 167 MeV. The given M L ,o and T* values correspond to the largest A point 
of each set. The discretization is characterized by N T = 512 and N s = 3 x 2 10 except for the points of set (a) at T — for 
A/Ap > 0.04 where N T — 3 x 512 was used. The step AL/T* was 5 for the cases (a) and (b) and 1 for case (c). 



Landau pole is high enough for the result not to de- 
pend too much on the value of the cutoff A. In order 
to study the sensitivity of the results on A we monitored 
the cutoff dependence of the relative change AQ(A) = 
(Q(A + AA) /(2(A) - 1) of a given physical quantity Q. If 
this quantity is regarded as a sequence for discrete val- 
ues of A, then in the ideal case where the convergence 
occurs the relative change not only tells us how sensitive 
is the physical quantity on the cutoff at some value A, 
but also how close it is to the convergent value. This is 
because if at some cutoff A the value |AQ(A)| < 10 _n , 
then one can say that Q(A) is within 10 _n+2 % from the 
converged value of the physical quantity Q. The problem 
is, of course, that strictly speaking the convergence would 
occur as A — > oo, but generally one cannot go above the 
Landau pole. Therefore, what is of practical relevance 
is whether the Q shows a plateau as a function of A be- 
low the scale of the Landau pole. We investigate this in 
Figure [5] for several parameter sets using the quantities 
</> and M 2 at different temperatures (the relative change 
is shown in percentage). One can see that we are clos- 
est to a plateau if the scale of the Landau pole is high 
and the temperature is low. The variation of the relative 
change with the cutoff shows that even when the scale 
of the Landau pole is approximately seven times larger 
than M^o, for practical purposes the result can be con- 
sidered compatible with a cutoff independent result, at 
least for temperatures not too large with respect to T pc . 



This result should however be interpreted with a pinch of 
salt since the fact that the plateau observed in Figure [5] 
extends up to the Landau scale is related to the fact 
that the physical quantities do not diverge at this scale, 
only the bare couplings. In higher order approximations 
where, due to a negatively quadratic growth of the self- 
energy at large frequency /momentum or to vertex- type 
resummations, one expects physical quantities to diverge 
at A p , it is less probable that a plateau can appear if the 
Landau scale is too lowF^I 

C. Phase transition 

In the chiral limit, in both the two-loop and the hy- 
brid approximations, the model undergoes a second order 
phase transition for those parameters of the (toJ/T! 2 , A*) 
plane, which are located above the T c = line of Figure[T] 
This is illustrated within the two-loop approximation in 
Figure® where we show the temperature evolution of the 
field expectation value, curvature masses and gap masses 
at the lowest available momentum. The inset shows that 



Some explicit calculation using an educated example reveals that, 
if A p is large enough, a plateau develops and even extends up to 
the very vicinity of A p . However, as A p is decreased, the plateau 
fades away. 
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FIG. 6. Illustration of the second order nature of the phase 
transition in the chiral limit within the two-loop approxima- 
tion. The inset shows the convergence of the static criti- 
cal exponents to the mean-field values, as the upper bound- 
ary of the fitting range is decreased. The parameters are: 
ml/T* = 0.124, A* = 22.277, while the discretization is char- 
acterized by A/T* = 100, N s = 3 x 2 10 and N T = 512 for 
T > 77 MeV and N T = 2 x 2 10 for T < 77 MeV. We also 
defined M L/T = M L/T (0, Aft). 



the three numerically determined critical exponents are 
compatible with the values /3 = 1/2, 7 = 1 and 5 = 3: 
thus the critical behavior of the system is characterized 
by mean-field-type critical exponents at this level of ap- 
proximation. This is expected, since they were already 
found to be of the mean-field type in the one-component 
case in Ref. @. 

One also sees in Figure [HI that Mt fulfills the re- 
quirement of Goldstone's theorem discussed in Sec. Ill Al 
around Eq. (fT7|> . as it vanishes in the broken phase 
and becomes degenerate with Ml in the symmetric 
phase. However, as a result of the truncation of the 
2PI effective action, Goldstone's theorem is violated by 
Mt(K = 0) (approximated numerically by Mt(0, Aft)) 
which is rather large since at small temperatures it is 
larger than 4>/2. We note however that Mr(0,Afc) is 
the smallest scale among M T (0, Aft), M L (0, Aft), M L 
and 4> an d that the size (in MeV) of the violation of 
Goldstone's theorem is quite constant with the temper- 
ature. These observations give good hope that higher 
order corrections will reduce uniformly the violation 
of Goldstone's theorem by the transverse gap mass. 
The restoration of Goldstone's theorem is expected 
because M = M in the absence of approximations. 
Our results indicate that the restoration could happen 
uniformly with the temperature. At large temperature 
both the degenerate curvature and gap masses increase, 
but a gap remains between them. This reflects the 
fact that the two-loop approximation is such that 
S^int/d^aS^bl^o ^ 2ST in t/5Gab\4>=o, where T int [cj), G] 



FIG. 7. The temperature evolution of the order parameter, 
the curvature and gap masses in the two-loop approximation 
where M L/ / T = M L / T (0, Aft). For <j> and Ml we also show 
for comparison the curves obtained in the hybrid approxima- 
tion, in which case we show the critical value of the field 
below which, at a given temperature, the effective poten- 
tial is not accessible. The parameters are: /T* = 0.04, 
A* = 17.39, h/T+ = 0.6. The discretization used in the two- 
loop case is characterized by A/T* = 55, N s — 3 x 2 10 and 
N T was increased for decreasing temperature from 512 used 
at T > 40 MeV to 2 x 2 10 for T € [25, 40] MeV and 4 x 2 10 
for T < 25 MeV. 



contains all the 2PI graphs contracted with the vertices 
of the shifted action S[tp —> tp+<f>] (see Ref. [l| for details). 

In the physical case (h ^ 0) the thermal transition is of 
an analytic crossover type. The temperature evolution of 
the order parameter is presented in Figure [7] for a set of 
parameters at which Ml,o ~ 360 MeV. The solid line is 
obtained in the two-loop approximation, while the barely 
distinguishable dashed line is obtained in the hybrid ap- 
proximation at the same values of the parameters. The 
inset shows that in the hybrid approximation 4>{T) is not 
a monotonous function of the temperature, for it shows 
a maximum at some value of the temperature. This re- 
flects an inconsistency of the hybrid approximation be- 
cause, as one sees in Figure [51 in the temperature range 
where 4>(T) > <f>o, the pressure is negative. In Figure 
the difference between the two approximations is more 
visible on the longitudinal curvature mass (the trans- 
verse ones differ very little because M| = h/cf>). The 
restoration of symmetry at high temperature is reflected 
by both the curvature and lowest momentum gap masses, 
as the corresponding longitudinal and transverse compo- 
nents approach each other. As in the chiral case, at large 
temperature there remains a gap between the curvature 
and gap masses. We note also the important difference 
between the values of M L (0, Aft) and M L > M L (0, Aft) 
at T = 0. It is clearly more convenient to use the cur- 
vature masses to parametrize the model since they allow 
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to reach higher sigma masses for the same values of the 
parameters. 



D. Thermodynamics 

We turn now to the study of the thermodynamic prop- 
erties of the model. To this end we compute the pres- 
sure by subtracting the value of effective potential at 
the minimum obtained at a given temperature from the 
value determined at zero temperature as described below 
Eq. (81) of Ref. @. The entropy density s = dp/dT is de- 
termined from the pressure through a numerical deriva- 
tive, while the energy is calculated as e = — p + Ts. 
Usually these quantities arc divided with appropriate 
powers of the temperature, but we choose to normal- 
ize them to the corresponding quantity calculated for an 
ideal gas of massless particles. As it was the case for 
N = 1 these three quantities when rescaled with the 
corresponding Stefan-Boltzmann limit agree with each 
other at that temperature where the interaction measure 
A = Td(p/T 4 )/dT = (e- 3p)/T vanishes. As discussed 
in Ref. @, this feature follows directly from the equa- 
tions, and can be seen in Figure [8j where we compare 
the dependence of these quantities on the temperature 
in the two-loop and hybrid approximations. The curves 
obtained in the two cases are indistinguishable above the 
pscudocritical temperature. At small temperature, how- 
ever, there are visible differences, and more importantly 
one can clearly see that the hybrid approximation is not 
consistent from a thermodynamic point of view, since at 
small temperatures it leads to negative pressure, entropy 
and energy densities. As already mentioned the temper- 
ature region where the pressure is negative is correlated 
to that where </>(T) > 0o- The inconsistency of the hy- 
brid approximation is displayed also by the heat capac- 
ity C — Td 2 p/dT 2 , which becomes negative for small 
temperature and by the square of the speed of sound 
c 2 = dp/de = s/C which has a singularity at the tem- 
perature for which C vanishes (see Figure [3]). Note, that 
for the two-loop case the temperature variation of c 2 is 
much milder than the one shown in Figure 8 of Ref. [32[ 
which was obtained in a Hartree approximation which 
included only the thermal effects and neglected the vac- 
uum ones. This is probably related to the fact that in our 
case Ml(T = 0) « 225 McV is smaller than the smallest 
value of the corresponding mass used there. 

It is visible in Figure [5] that at high temperature the 
pressure normalized to the Stefan-Boltzmann limit de- 
creases with the temperature. This is the consequence 
of the fact that, as one can see in Figure at high T 
the masses of the excitation grow linearly with T and 
therefore a high temperature expansion is less and less 
accurate with increasing T. 

In the two-loop approximation we also tested the de- 
pendence of the pressure on the renormalization scale 
T+, by choosing two points in the parameter space 
which belong to a line of constant physics, along which 
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FIG. 8. The temperature dependence of the pressure, energy 
density and entropy density normalized to the corresponding 
Stefan-Boltzmann limit obtained for the same parameter set 
of Figure [7J in the the two- loop (continuous lines) and hybrid 
approximations (dashed lines). The inset shows the variation 
of the interaction measure with the temperature. 
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FIG. 9. The square of the speed of sound Cg and the heat 
capacity C obtained for the same parameter set of Figure [7J 
in the the two-loop (continuous lines) and hybrid approxi- 
mations (dashed lines). The left and right axes of the plot 
correspond to and C, respectively. 



M L / Tj o(0, Afc) and </>o are constant, and for which the 
difference in Mt,o(0, Afc) was maximal, that is around 
10%. The difference in the value of T* corresponding 
to these two points was around 10% and although the 
difference of Ml,o(0, Afc) was around 30%, the maximal 
difference in the pressure was around 10% and was ob- 
served at temperatures smaller than T pc . 
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VI. CONCLUSIONS 

We studied numerically the thermal phase transition 
of the renormalized O(N) model, both in a genuine 3>- 
derivable approximation in which the effective action is 
truncated at two-loop level and in a hybrid approxima- 
tion in which the effective potential and the field equation 
derived from it are evaluated with a lower level, Hartrcc- 
Fock-type transverse (pion) and longitudinal (sigma) 
propagators. In the first case the self-consistent propaga- 
tor equations were solved iteratively in Euclidean space 
using 3D cutoff regularization and the method of Ref. Q , 
which by a combination of adaptive numerical integra- 
tion and fast Fourier transforms ensures a very accurate 
evaluation of the convolution-type integrals. In the hy- 
brid approximation one obtains explicitly finite equations 
which are much simpler to solve. 

In the chiral limit the phase transition turns out to 
be of second order in both approximations studied. On 
the one hand, this means that the higher level truncation 
considered in this work represents an improvement over 
the Hartree-Fock approximation which is known to yield 
a first order phase transition in the chiral limit. On the 
other hand, we have a clear indication that the impor- 
tant improvement over the Hartree-Fock level occurs in 
the field equation and is related to the inclusion of the 
setting-sun diagram. In the case of an explicit break- 
ing of the chiral symmetry the transition is an analytic 
crossover. 

As long as one is interested in the temperature evolu- 
tion of the expectation value of the field, curvature and 
gap masses the hybrid approximation can be regarded 
as a good approximation of the two-loop ^-derivable ap- 
proximation. In the chiral limit this is not true for the 
entire parameter space, as one has to restrict its appli- 
cation to those parameters where the longitudinal curva- 
ture mass does not change abruptly with the parameters. 
However, the thermodynamic study revealed its inconsis- 
tency at small temperatures for it leads to negative pres- 
sure, entropy density and energy density. In fact, this 
feature is also related to the nonmonotonic behavior of 
the field expectation value at small temperature, where 
it first increases with increasing temperature. 

We have seen that for N = 4 it is possible to achieve 
a realistic parametrization of the model, in which the 
zero temperature sigma mass, obtained as the longitu- 
dinal eigenmode of the curvature tensor, could be fixed 
to values around 460 MeV, while keeping the scale of 
the Landau pole at around 3.4 GeV. This scale is large 
enough for the results to be considered practically inde- 
pendent of the cutoff used, at least for the approxima- 
tion considered here and for temperatures not too large 
with respect to the crossover temperature. The values of 
the sigma mass which can be obtained within the two- 
loop 2PI approximation are larger than those found in 
the next-to-leadin g or der of the 1/N expansion in the 
1PI formalism [28U29j). and the scale of the Landau pole 
proved also larger. However, in the approximations stud- 



ied here, there is a significant difference between the cur- 
vature masses and the gap masses. It is expected that in 
approximations where the effective action is truncated at 
higher orders this discrepancy will diminish and then the 
question is raised whether this will affect the maximum 
value of the sigma mass achievable from the curvature 
mass. In this respect, it will be interesting to investigate 
whether the possibility of a realistic parametrization of 
the 0(4) model persists in the 2PI formalism at higher 
order truncation levels and also what will be the case in 
the linear sigma model with three flavors at two-loop and 
higher truncations levels. Also, in higher order approxi- 
mation, especially in those involving vertex-type resum- 
mation, as the 2PI-1/JV expansion, the presence of the 
Landau pole is a more severe problem, because it influ- 
ences the renormalized quantities due to the divergence 
of the vertex function. In this case the cutoff insensitiv- 
ity needs a careful reexamination and it is more probable 
that the scale of the Landau pole has to be kept fur- 
ther away from the physical scales than in the two-loop 
approximation discussed here. 
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Appendix A: Four-point functions 

As it is discussed at length in Ref. [l[ there exist three 
different definitions of the four-point function which do 
not match exactly in a given truncation. Here we shall 
consider these distinct definitions at <j) = in which case 
we can use G^~° = 5abG<p = o- One possible definition i j^l 

-0-0 -0-0 1 /* T -0-O "2 "0-0 

V a b,cd = ^ab,cd ~ 2 J Q ^aM«^0wK»,a!' (^-1) 

with 

A = o = 4<5 2 $ fA2l 

ab ' cd ~ SGab(K)SGcd{Q) 0=o ' y ' 

The kernel A^° d has the structure Af^° d = A^ Q S ab S cd + 
AgjfoAd + S ad S bc ) with A^ B) = X^ B) /3N. The 



In the approximation at hand, the kernel, and in turn the four- 
point function, are momentum independent, hence the absence 
of momenta in our notation for the kernel and the four-point 
function. 
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four-point function V^ b c ° admits the same decomposition 
and its components obey the coupled set of equations 



fM) _ a (A) —I \^ A > f' 1 (n\v w 
v <f,=a ~ A 0=o _ ~K / Vo°*=oiyJVo 



N 



(A) r,2 



r{A) 



(-1) 



V <^.=0 — Ji 0=O 



^fJ GU(Q)V^. (A4) 



(B) 



By expanding this system perturbatively, it is pretty ob- 
vious that and are not equal and thus that 
V a b,cd is not crossing symmetric. In particular, the diver- 
gent part is not crossing symmetric, which explains the 

need for two distinct bare couplings A ^ and X^ B \ In or- 



der to solve the system of equations (|A3j) - ()A4j) . it is con 



venient to consider the combinations A 



2A<*> = X { n NA+2B) 



is then easily proven that 



(CO 

0=0 



(CO 
<f>=o 



NA 



NV„ 



(A) 



(A) 
0=0 



0=0 



(C) 
=0 



A 



(CO 
>=0 



1 



A£ ) G^ (Q)^, (A5) 



(C) 



which shows that the combination V^_q diagonalizes the 
system (|A3j) - (|A4|) . This diagonalization is in one-to- 
one correspondence with the one we used in the case 
of the Hartree gap equations in Sec. IIIII In fact the 
combinations of bare couplings A^ = A /(3AT) and 



XK"^^' /(3iV) are precisely those which 



J^(NA+2B) _ ^(NA+2B) 

appeared in the diagonalized form of the Hartree-Fock 
gap equations. We obtain 



1 



V 



(B) 
=0 



3N 

7(B) 



■B[^=o](0) 



37V 



V, 



(C) 
=0 



(NA+2B) 



X 



:^[G0= O ](O). 



(A6) 
(A7) 



which we use in the main text to derive the expressions 
for A and Xq NA+2B ^ from the renormalization and con- 
sistency conditions. A second definition of the four-point 
function which depends on one momentuirF^l is 



y£2m - k:u k ) § / K::LG'U(Q)v:=: d (Q) 

z Jq 

= - \ f v£: v Gl= (Q)At:° d (Q) , 



(A8) 



with 



25 3 $ 



= 



8G ab {K)8(j) c 8(l) d 



k<£2 Q 5abScd + A 0=o y>ac8bd + 5 ad 5bc 



(A9) 



(A3) and 



A 



A 



(A) 

= 

(B) 
1=0 



1 

3N 
1 
3N 



X^- 

X (B) 



3N 
W + 6 
6N 



XlBiG^oKK) 



(A10) 



A^[G0 =o ](K) . (All) 



Once again the appropriate combination of components 



NV 



=0 -r 2V^ = g leads to a system of decoupled equations 

( B) (C) 

for V^ = g and V^_q which is suited in particular to extract 

the expressions for the bare couplings X^' B \ According 
to Ref . [lj , the third possible definition of the four-point 
function is given by 

<5 4 7 



V. 



i=0 
abed 



1=0 



__1 

3N 

1 

~ 2 
1 

~ 2 
1 

~ 2 



A^(0)G^ (Q)l^(Q) 



aL:L(Q) g 0=o(Q)CX(Q) 



0=0, 



(A12) 

Using the tensor decomposition of V^=o, we check that 
VabZd = ^0=o(^a6<5cd + SacSbd + SadSbc) (in other words 



^0=o = ^0=o an( I tnus ^0=o nas the crossing symmetry) 
with 

V0=o 



r(B) 



X4_ 

3N 



1 

2N 



^ (Q)GU(Q)v^>(Q) 



2 1 



1 

A 



A 



S(O)^ (Q)^(Q), 



(B), 



(A13) 

where we note that the contributions from the B and C 
components of A^ =0 an d V^—o contribute independently 
to 1^0=0 ■ This expression is used in the main text in order 
to obtain A4. 



Appendix B: Hybrid extras 



18 In principle, a four-point function depends on three independent 
momenta. The four-point functions we consider here are taken 
for particular values or configurations of their external momenta 
and can thus depend on fewer variables. 



In this section we give some expressions encountered 
in the hybrid approximation and discuss their particular 
aspects in some details. 
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1. Expression of V^=o 



\V {A) -A iA) ] ~ 2\V (B) - A (B) 1 
l v <h=o y V=oJ z L^<*=o 7 V=oJ 



(B2) 



As mentioned in the main text, the four-point function 
V0=o is modified in the hybrid case. Following the same 
strategy as in the previous section, we arrive at 



y0 = O _ 

abed 37V V 

f 

~ 2 
f 
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Sabred + SacSbd + S a dSb< 
T 



^td^uv (Q)G^, = q {Q)V UVi bc 

(Afrl(Q)-Af=°J^ =0 (Q)Kt 



(A2=S„(Q) - At:l)G| =0 (Q)t^t d 



(Aad,L(Q) ~ A ad ° !J )G| =0 (Q)^ i6 c . 

(Bl) 



Once again V^ b J = V^ =0 (5 o 6^cd + <5ac<5&d + <5 ad <5 6c ), with 
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3N 
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>=oJ 



from which we deduce Eq. (|40 



2. Finiteness of C and C 

Let us show that C and C defined in ([51]) and (f54"|) are 
finite. One possibility is to compute the divergent part 
of the setting-sun sum integrals and check that C and C 
are free of divergences. As far as C is concerned this was 
done within dimensional regularization in Appendix B of 
Ref. 0. We can treat C along similar lines. With the 
exception of the setting-sun integral <S[Gl; Gt; Gt], all 
the integrals needed to obtain the finite expression of C 
are given in Ref. [§] . Using the method of Ref. [33| , one 
can obtain for this setting-sun integral a decomposition 
in terms of zero, one and two statistical factors analo- 
gous to Eq. (B5) of Ref. @- From that point on, the 
calculation of 5[Gl; Gt; Gt] an d C[Gl, Gt, G*] parallels 
that of S[G] and C[G, G+] performed there and uses the 
vacuum part of the setting-sun integral with two differ- 
ent masses. For the part with no statistical factors one 
has to expand the factors of the product 7f[G]-B*[G*](0) 
to 0(e) because both contain 1/e divergences. Using for 
iS*-°) the expression given in Sec. 3 of [34j one obtains 
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where ^i(x) = c? 2 r(a;)/(ix 2 is the trigamma function, z = Ml/(4M|), and the function 3?(z) is defined as 
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with Ch(x) = - J d6 ln(2 sin (9/2)) being the Clausen function and £(z) = y/l — l/z. Note that lim $(z) 
o 2-5-1 



The part with one statistical factor reads 
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7T 

7! 



M 2 



if z > 1, 
if z < 1, 



Sin 2. 



(B5) 
(B6) 

(B7) 
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with Q — yfl — Finally, the part with two statistical factors is 



C^[Gl,G t ,G,} = 2T^[G T ]B[ 1 \G,}(0) - 1 3^ 2 )[G l ;Gt;G t ] -S (2) [G l ] - 2<sf 2) [G,] - 2(M 2 - m 2 ) d5 f 

3 amf 



(B8) 



where 



32tt 4 



dp / dkpk 



n T (e T {k)) 
e T (k) 



n T (e T (p)) ln 44(fc)e 2 (p) - (M 2 - 2M 2 + 2fcp) 2 



:(P) 



44(fc)4(p) - (Af 2 - 2M 2 - 2fcp) 2 



nr(e L (p)) , 4 £ |(fc) £ 2 (p) - (M 2 - 2fcp) 2 
e L (p) 4 £ 2 (fc) £ 2 (p)-(A7 2 +2fcp) 2 



(B9) 



with n T (e) = l/(exp(e/T) - 1), £ T/L (A:) = fc 2 
and all other integrals are given in Ref. |{|. 



M 2 
T/L' 



We can also discuss the finiteness of C and C using 
another explicit method, which is less calculations! and 
shows clearly the different roles played by A21 and <5A2 n i- 
Let us illustrate it in the case of C . Using the results of 
Ref. [331 to write 



T{G}=Tt= [G] 



(BIO) 



and 



S[G] =S t=0 [G] + 3 f^Yl BffllQ') 
J i Z£ i <j=±i 

+ 3 /? 1 It" E G{Q„ + K T UB\\) 

with / = / d 3 q/(2Tr) 3 and where B[G](Q a ) denotes the 
analytical continuation of the bubble sum integral to real 
values of the frequency followed by its evaluation on shell 
Q = (iu) n ,q) -> Q a = (q = ae q + iO + ,q). Plugging this 
back into the first line of Eq. (|5Tj) , we obtain 



C[G,G,] = A^_ 



finite 



T=0 



<JG(Q)[B*[G*](O)-0*[G*](Q)] 



gi 2 [b4gj(o)-b[g](q c 

^ (T — ±1 

/S- E k[G,](0)-S,[G,](Q.) 
J ? ze ? CT =±i 



(B12) 



where 5G = G - G*, e g = vV + M 2 and e* = 
\Jq 2 + mf. The role of 5A2ni is thus to remove the subdi- 
vergences in such a way that the integrands in Eq. (|B12|) 
are finite. In addition the last two integrals are finite 
due to the presence of the thermal factors. This is not 
the case for the zero temperature integral because the 



factor 5G does not decrease fast enough in the UV. How- 
ever the term involving the ratio A21/A0 has precisely the 
same form as this integral. Separating the T = part of 
A21/A0 and combining it with the rest we obtain 



C[G,G*] = 



T=0 



G r (Q) [S*[G*](0)-B*[GJ(Q)] 
[b*[G*}(0) - B[G](Q a ) 



u=±i 



/fe E [B*[G*](0)-B*[G 



finite . 



(B13) 



where G r = SG + (M 2 - m^G 2 . = (M 2 - m 2 ) 2 GjG 
decreases fast enough to make the corresponding integral 
convergent. A similar proof can be given for C. 



3. Explicitly finite field equation 

We give here the explicitly finite form of the field equa- 
tion. Differentiating Eq. ((52")) with respect to (j> and tak- 
ing into account the implicit dependence on <j> of the gap 
masses one obtains 



Ml 



3N ( 



A 



18iV 2 



Cn[Gl, Gt, G*] 



+4? ((JV + 8)V[G L , GJ + (JV - 1)T> L [G L , G T ] 



dA^ 



+<^ 2 (^-1)Pt[G l ,Gt,G,] 



0, 



(B14) 



where, in order to give a compact expression, we have 
already used in the first two terms, coming from the 
Hartree-Fock part, the solution of the two linear equa- 
tions for the derivatives of the gap masses which are ob- 
tained from Eqs. ((36} and ((47]) as 



dM 2 



d-b 
ad — cb 



dM\ 
d4> 2 



ad — cb 



(B15) 
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with c = B[G h ](0) - B*[G*](0), a = c + 2JV/A*, 6 = 
(iV-1) [B[Gt](0)-B*[G*](0)]/3, and d = 36(iV+l)/(iV- 
1)+6AT/A*. Both 6 and c are finite because the divergence 
of the perturbative bubble integral does not depend on 
the mass. We note that for N = 1 one has 6 = and 
dMl/d4> 2 = 1/a 
the field equation simplifies to 



therefore with the notation M 2 = M 2 



= 



M 2 — — 



A 2 



yC[G,G. 
V[G, G*] 



A* x + [B[G](0)-£*[G*](0)]/2 



(B16) 



In Eq. pli]) . C N is defined in Eq. ((SSJ) and X>[G L , GJ 
whose explicit expression is given in Eq. (B12) of Ref. Q 
is obtained from C[Gl,G*], defined in (fBTj) . upon diffcr 



cntiation with respect to M^. Finally, 2?l[Gl,Gt] and 
2>t[Gl, Gt, G*] are obtained from C [Gl , Gt , G*] given 
in Eq. ((54)) upon differentiation with respect to M 2 and 
M|, respectively. 



-l/e 



Wo(x) 




-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 

FIG. 10. The two real branches of the Lambert W func- 
tion. The upper branch Wo (a;) (dashed) is defined for x £ 
[— l/e, oo) and the lower oneW-i(i) (solid) for £ € [— l/e, 0). 



In order to determine the 4>c,o curve in the (m 2 /T 2 , A*) 
plane, discussed below Eq. (f73")) , we have to take the limit 
T -> and Af T -> in Eq. (|Bl4]) . This is straightfor- 
ward, except for the last term where 2?t, which at T = 
comes only from C*- -*, diverges and dM^/dcj) 2 vanishes in 
this limit. Working out the limit one obtains 



lim 



~D t[Gl, Gt, G*] 



dM 2 " 1 



i 



hl 



i\/f 



dcj) 2 

,0 1 



T=0 
(1) 



26i 1J [G,](0) 



(JV+2)A 
3JV 



i 

16tt ; 



M 2 



AT- 1 



.(B17) 



As a last remark, note that there is no equivalent of 
the 4> c fi line in the N = 1 case, since this requires </>o > 
and Mq = 0, and by examining Eq. (|B16|) at T = 
one sees that the two conditions cannot be satisfied 
simultaneously This is because, for Mq — > 0, C is finite, 
but T> diverges as ln 2 (Mo/™*), while the denominator 
diverges only as ln(Mo/m*), meaning that in this limit 
(j)Q = is the only solution of the field equation. 



4. Solution of the gap equation at (f> — and T = 



We show here that the solution(s) of Eq. (|42|) can be 
given at T = in terms of the two real branches of the 
Lambert W function defined to be the multivalued in- 
verse of the function w i— > iiie 1 " = z, for m; complex. This 
function verifies W(z) exp(W(z)) = z for any complex 
z. The real branches of the Lambert W function are de- 
picted in Figure 1101 The upper branch is usually called 
Wo (a;) and the lower one W_i(x). 

At T = one can rewrite Eq. (g2J) as 



M^ =0 T=0 In 



M 2 

o 



-— , (B18) 
a* 



where a* = (N + 2)\+/(96ir 2 N) 1 = -1 + (N + 

2)Aji3i 1) [G*](0) - 1/(16tt 2 )]/(6A0 and G* is defined in 
Eq. (|67p . With a few algebraic manipulations (exponenti- 
ation and multiplication by — G*/ (a*M? =0 T=0 )) and us- 
ing the definition of the Lambert function, one expresses 
the solution of Eq. (|B18|) as 



M. 



0=0, T=0 



W 



(B19) 



For C* > (points below the T c = line of Figure □]) 
the argument of W is negative and one sees by looking 
at Figure [10] that for M^, = o.t=o one has no solution if 



c. 



< — l/e and two solutions if 



< 



— l/e, one smaller and one bigger than M e = 2Ap St /e, 
where A p st is the accurate estimation of the Landau pole 

given in (|68|) (the two solutions merge when C ir e b *^ a * = 
m+a+fe). The lower scale solution is given by the lower 
branch W-i of the Lambert function[3 The larger scale 
solution is given by the upper branch Wo- For G* < 
one has one solution, bigger than M e and given by Wo. 

In conclusion, to obtain the T c = curve, defined as 
M^,=o,r c =o = 0, we have to take the solution (|B19|) given 
by the lower branch of the Lambert function and use it 
in Eq. (|70p. which can be solved only numerically. For 
small negative arguments W_i can be evaluated using 
the asymptotic series given in Ref. [Hj]. We mention 
finally that the upper branch plays a role because we 



19 We know that at T c = one has M^_q j,^_ = 0, and this can 
be obtained only with W_i, which diverges negatively when its 
argument vanishes. The use of the other branch would give a 
finite value, because ~V\!q(x) = x + 0(x 2 ), for small x. 
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have considered the rcnormalized gap equation in its 
continuum limit. If we would consider it in the presence 
of a finite 3D cutoff, the solutions corresponding to the 
upper branch would be absent, see Rcf. ha. 



Appendix C: Tensor decomposition 

Let us first consider a real symmetric tensor U% b = U^ a 
such that for any rotation R £ SO(N) 

U^ = R ac R bd ut d . (CI) 

The case <f> — is easily treated. For any AgIwc have 

([/f=° - XS ab )R bc = Rab(Ut° - A<M . (C2) 

Since the fundamental representation of SO(N) is irre- 
ducible, it follows from Schur's lemma that U — Al is 
either or invertible. If we chose A to be an eigenvalue 
of U^ =0 (there exists at least one eigenvalue since U^ =0 
is real and symmetric), {7* =0 — Al cannot be invertible 
and thus [/^ =0 = Al. For reasons that will appear below, 
the case <j> ^ requires that we distinguish N = 2 from 
N > 2. Let us consider the case N > 2 first. Since is 
real and symmetric, it is diagonalizablc, that is it admits 
N linearly independent eigenvectors. Let us consider an 



eigenvector which is not collinear to <f> (such an eigen- 
vector exists for N > 2). If R is a rotation that leaves 4> 
invariant, we have from (|CI[) : 

U&Rbcut = RabUtut = X^Rab4 , (C3) 

which shows that we have indeed at least N — 1 linearly 
independent eigenvectors corresponding to the eigenvalue 
X 4 . If the remaining eigenvector corresponds also to A*, 
we have U 4, = XH with X R<I > = A* that is with A* a func- 
tion of 4> 2 only. If the remaining eigenvector corresponds 
to another eigenvalue ^ X 4 , it has to be collinear to 
<f>. If it were not, we could construct N — 1 linearly in- 
dependent eigenvectors, different from the previous ones 
and it would follow that TV > 2N — 2 that is TV < 2, 
which contradicts our assumption N > 2. From this and 
the fact that the subspaces are orthogonal to each other, 
it follows that 

U* ab = + X<t>Pl, (C4) 

with A^ and functions of cf> 2 only. The case N = 2 is 
particular because the solution to Eq. (|Cip is much more 
general than (|C4[) . In fact, because SO(2) is Abelian, 
any tensor of the form 

U? lb = R ac R bd (//P £ L d + \+P%) , (C5) 

with R £ 5*0(2), obeys Eq. (|C1|> . The converse can also 
be proven to be true. By imposing that Eq. (|C1[) holds 
not only for any R £ SO (2) but also for any R £ 0(2), 
one recovers the form (|C4[) . 
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